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ABSTRACT 

The afterglow emission from gamma-ray bursts (GRBs) is usually interpreted as synchrotron radia- 
tion from electrons accelerated at the GRB external shock, that propagates with relativistic velocities 
into the magnetized interstellar medium. By means of multi-dimensional particle-in-cell simulations, 
we investigate the acceleration performance of weakly magnetized relativistic shocks, in the magne- 
tization range < a < 10 _1 . The pre-shock magnetic field is orthogonal to the flow, as generically 
expected for relativistic shocks. We find that relativistic perpendicular shocks propagating in electron- 
positron plasmas are efficient particle accelerators if the magnetization is a < 10~ 3 . For electron-ion 
plasmas, the transition to efficient acceleration occurs for a < 3 x 10 -5 . Here, the acceleration process 
proceeds similarly for the two species, since the electrons enter the shock nearly in equipartition with 
the ions, as a result of strong pre-heating in the self-generated upstream turbulence. In both electron- 
positron and electron-ion shocks, we find that the maximum energy of the accelerated particles scales 
in time as e max oc i 1 / 2 . This scaling is shallower than the so-called (and commonly assumed) Bohm 
limit Emax oc t, and it naturally results from the small-scale nature of the Weibel turbulence generated 
in the shock layer. In magnetized plasmas, the energy of the accelerated particles increases until it 
reaches a saturation value e sa t/"/o m iC 2 ~ c -1 / 4 , where jorriiC 2 is the mean energy per particle in 
the upstream bulk flow. Further energization is prevented by the fact that the self-generated turbu- 
lence is confined within a finite region of thickness cx c -1 / 2 around the shock. Our results can provide 
physically-grounded inputs for models of non-thermal emission from a variety of astrophysical sources, 
with particular relevance to GRB afterglows. 

Subject headings: acceleration of particles - cosmic rays - gamma-ray burst: general - pulsars: general 
- radiation mechanisms: non thermal - shock waves 



1. INTRODUCTION 

The external shocks of gamma-ray bursts (GRBs) are 
often invoked as effi cient sites of acc eleration for protons 
and electrons 



(e.g. 



Waxman 2006). Shock-acceleration 



of protons might explain the flux of Ultra High Energy 
Cosmic Rays (UHECRs) observed with energie s in excess 
of 10 20 eV by the Pierre Auger Observatory (Abraham 



20101. Synchrotron emission from the shock-accelerated 
electrons powers the GRB afterglow emission, which is 
usually detected in the X-ray, optical and sometimes 
radio bands, and recent ly up to sub-GeV energies by 
the Fermi telescope (e.g., Ackermann||2010 De Pasquale 
2010) IGhisellini et al.||20ltf 



The external shocks in GRB afterglows are believed to 
be relativistic shocks propagating in weakly magnetized 
electron-proton plasmas, either the interstellar medium 
(ISM) or the wind of the progenitor star. If the ISM 
number density is n = uq cm' 3 and the magnetic field is 
Bism = 3 -Bism,-5.5 A*G, the ISM magnetization will be 
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where m.j is the proton mass and c is the speed of light. 
For relativistic shocks, the mean field in the post-shock 
frame will be mostly transverse to the flow, due to shock 
compression and to the effect of the Lorentz transfor- 
mations. For perpendicular shocks (i.e., with the field 



orthogonal to the flow), the magnetization in eq. is a 
Lorentz invariant, and independent of the shock radius 
(provided that n and -Bism are constant in radius) rj 

The acceleration process at the external shocks of 
GRBs is thought to be governed by the Fermi mech- 
anism, where particles stochastically diffuse back and 
forth across the shock front and gain energy by scattering 
from mag netic turbulence embedded in the converging 



flows 



e^ 



"Blandford fc Ostrike r | [l978l [B^fl|[l978l |Drury 
19831 |Blandford fc Eichler||l9¥7f .^TrIe highly nonlinear 
physics of the Fermi process - where the magnetic turbu- 
lence that mediates the particle acceleration is generated 
by the particles themselves - can only be addressed self- 
consistently by means of first principle parti cle-in-cell 
(PIC) simu l ations. By using PIC simulations, Sironi & 
Spitkovsky (2009a) (in electron-positron flows; hereafter, 
SS09) and |Sironi fc Spitkovsky (2011b I (in electron-ion 
flows; hereafter, SS1I) have demonstrated that parti- 
cle acceleration is suppressed in a ~ 0.1 perpendicular 
shocks. Here, the self-generated turbulence is not strong 
enough to permit efficient injection into the Fermi pro- 
cess. On the other hand, in the extreme case of un- 
magnetized shocks (i.e., a = 0), PIC simulations have 
shown that a non-thermal tail of shock-accelerated par- 
ticles is generated, as a self-consistent by-product of the 

1 The independence of a from the shock radius also holds for a 
wind profile of the external density, under the assumption that the 
field is primarily toroidal. 
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Fig. 1. — Simulation geometry. For our 2D simulations, the com- 
putational domain is in the xy plane, with periodic boundary con- 
ditions in the y direction. The incoming flow propagates along 
—x, and the shock moves away from the reflecting wall (located 
at x = 0) toward +x. The magnetic field carried by the upstream 
flow (red arrows) is perpendicular to the simulation plane. 



shock evolution jSpitkovsky|2008a|b||Martins et al.|2009 
Haugb0lle|[Ml] ) 



in this work, we investigate the physics of weakly mag- 
netized perpendicular shocks for magnetizations in the 
range < a < 10 _1 . We explore how the accelera- 
tion efficiency, namely the fraction of particles and en- 
ergy stored in the non-thermal tail, depends on the flow 
magnetization, in electron-positron (f|3) and electron-ion 
(;fll) flows. We show that shocks propagating in electron- 
positron plasmas are efficient particle accelerators, if the 
magnetization is a < 10 -3 . A smaller threshold is found 
in the case of electron-ion plasmas (a < 3 x 10~ 5 ). Here, 
the acceleration process proceeds similarly for the two 
species, since the electrons enter the shock nearly in 
equipartition with the ions, as a result of strong pre- 
shock heating in the self-generated turbulence. 

In all the cases where the Fermi process is efficient, we 
follow the temporal evolution of the upper cutoff of the 
non-thermal tail of accelerated particles, finding that the 
maximum energy scales as E max oc i 1 / 2 , for both electron- 
positron and electron- ion flows. This is in contrast with 
the ad hoc prescription given by the so-called Bohm scal- 
ing (namely, e max cx t), which is often employed in the 
literature, for lack of a better choice. In magnetized plas- 
mas, we find that the energy of the accelerated particles 
increases up to e sat / / 7om i c 2 ~ o --1 ^ 4 , where 7 m i c 2 is 
the mean energy per particle in the upstream bulk flow. 
Further energization is prevented by the fact that the 
self-generated turbulence is confined within a finite re- 
gion of thickness oc cr -1 / 2 around the shock. 

With our PIC simulations, we are then able to provide 
a physically-grounded scaling for the acceleration rate 
in weakly magnetized relativistic perpendicular shocks. 
Our results can be easily incorporated in models of non- 
thermal emission from astrophysical sources, as we illus- 
trate in f|5] We finally summarize our findings in [J6j 

2. SIMULATION SETUP 

We use the three-dime nsional (3D) elec tromagnetic 
PIC code TRISTAN-MP dSpitkovsky|j2005| , which is a 
paral lel version of the public code TRISTAN ( Buneman 
|1993| that was optimized for studying relativistic coili- 
sioniess shocks. Our simulation setup parallels SS09 and 
SS11 very closely, which we repeat here for completeness. 

The shock is set up by reflecting a cold "upstream" 
flow from a conducting wall located at x = (Fig. [IJ. 
The interaction between the incoming beam (that prop- 
agates along — x) and the reflected beam triggers the 



formation of a shock, which moves away from the wall 
along +x. This setup is equivalent to the head-on col- 
lision of two identical plasma shells, which would form 
a forward and reverse shock and a contact discontinu- 
ity. Here, we follow only one of these shocks, and replace 
the contact discontinuity with the conducting wall. The 
simulation is performed in the "wall" frame, where the 
"downstream" plasma behind the shock is at rest. 

We perform simulations in both 2D and 3D compu- 
tational domains, and we find that most of the shock 
physics is well captured by 2D simulations. Therefore, 
to follow the shock evolution for longer times with fixed 
computational resources, we mainly utilize 2D runs, but 
we explicitly show that our 2D results are in excellent 
agreement with large 3D simulations. For both 2D and 
3D domains, all three components of particle velocities 
and electromagnetic fields are tracked. In our 2D sim- 
ulations, we use a rectangular simulation box in the xy 
plane, with periodic boundary conditions in the y di- 
rection (Fig. nl). In 3D, we employ periodic boundary 
conditions both in y and in z. Each computational cell 
is initialized with four particles (two per species) in 2D 
and with two particles (one per species) in 3D. We have 
performed limited experiments with a larger number of 
particles per cell (up to 8 per species in 2D), obtaining 
essentially the same results. 

We investigate the physics of both electron-positron 
and electron-ion shocks. In both cases, the relativistic 
skin depth for the incoming electrons (c/w pc ) is resolved 
with 8 computational cells, and the simulation timestep 
is At = 0.056 u;",, 1 . Here, w pc = (47re 2 n e /7 m e ) 1 / 2 is the 
relativistic plasma frequency for the upstream electrons, 
with number density n e (measured in the wall frame) and 
bulk Lorentz factor 70. For electron-positron flows, the 
plasma frequency of the upstream positrons is w p j = oj pc . 
In the case of electron-ion shocks, we typically employ a 
reduced mass ratio m i /m e = 25, which allows to follow 
the shock evolution for sufficiently long times (in units of 
the inverse ion plasma frequency uj~ 



^mi/m e Wpc 1 ), 
while still clearly separating the ion and electron dynam- 
ical scalesF] As we show in Appendix [A we obtain es- 
sentially the same results when using higher mass ratios 
(we have tried up to rrii/m e = 1600, approaching the re- 
alistic value mi/m e ~ 1836), which suggests that a mass 
ratio mj/m e = 25 is already "large" enough to capture 
the acceleration physics in our electron-ion shocks. 

For electron-positron shocks, our computational do- 
main in 2D is typically ~ 128c/w P i wide (corresponding 
to 1024 cells), and in 3D the transverse size of the box 
amounts to ~ 64c/w p i, or 512 cells. We have tried with 
2D boxes up to three times as wide, finding essentially 
the same results. In electron-ion shocks, we choose a box 
with a width of 1024 cells in 2D and 512 cells in 3D, which 
correspond respectively to ~ 26c/w p i and ~ 13c/w p i, for 
our reference mass ratio mi/m e = 25. When scaling to 
higher mass ratios, we choose a 2D computational do- 
main with 1024 transverse cells for m,/m e = 100, with 
2048 cells for m./rrie = 400, and with 4096 cells for 

2 We remark that in the following we use the same symbol ui p i to 
indicate the relativistic plasma frequency of positrons (in electron- 
positron flows) and of ions (in electron-ion shocks). Similarly, 
we indicate with m; the mass of the positively-charged particles 
(positrons or ions, depending on the flow composition). 
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mi/m e = 1600, so that in each case the width of the 
box amounts to ~ 13c/w p i. 

The incoming plasma is injected through a "moving 
injector," which recedes from the wall along +x at the 
speed of light. The simulation box is expanded in the 
x direction as the injector approaches the right bound- 
ary of the computational domain. This permits us to 
save memory and computing time, while following the 
evolution of all the upstream regions that are causally 
connected with the shock. The final time of our sim- 
ulations is chosen such that we can confidently predict 
the subsequent evolution of the shock. Our longest 3D 
simulation ran for 50,000 timesteps, corresponding to a 
box with 22, 500 cells along x, while in 2D we were able 
to evolve a computational domain with 3072 transverse 
cells up to 400, 000 timesteps, at which point the box 
extends over 180, 000 cells in the x direction. 

The cold incoming stream is injected along — x with 
bulk Lorentz factor 70 = 15. As we show in Sj3] and 

four results are nearly the same for higher values 
the Lorentz factor (we have explored from 70 = 3 
up to 70 = 240), modulo an overall shift in the en- 
ergy scale. The upstream flow is seeded with a back- 
ground magnetic field B , which we parameterize in 
terms of the ratio of magnetic to kinetic energy density 
a = Bq I ' AirjomiriiC 2 , where (=?v) is the number den- 
sity of incoming positrons or ionsrl We investigate the 
regime of weakly magnetized flows, with magnetizations 
lower than the a = 10 _1 case explored by SS09 and SS11, 
down to th e limit of unmagnetized shocks discussed by 
Spitkovsky ( 2008a|b ). The regime of magnetizations we 
study here may be relevant for e xternal shocks in GRB 
afterglows, as we discuss in §5 1 1 

The angle between the shock direction of propagation 
+x and the upstream magnetic field Bo is taken to be 
9 = 90° in all cases, i.e., we only focus on "perpendicu- 
lar" shocks (Fig. As discussed in fjl| this is the most 
relevant configuration for the relativistic shocks of GRB 
afterglows in the frame of the post-shock medium, due 
to Lorentz transformation effects and shock compression. 
We refer to SS09 and SS11 for a complete investigation of 
the dependence of the shock physics on the field obliquity 
9, from 9 — 0°, which corresponds to a parallel shock, 
with magnetic field aligned with the shock normal, up to 
9 = 90°, i.e., a perpendicular shock, with magnetic field 
along the shock front. In the upstream medium, we also 
initialize a motional electric field Eq = — /3 x Bo, where 
(3 = — /3o x is the three-velocity of the injected plasma. 
In our 2D experiments, we choose a magnetic field or- 
thogonal to the simulation plane, as shown in Fig. [l] By 
comparison with 3D simulations, we have verifiea that 
this is the field geometry most appropriate to capture 
the 3D physics of the shock, in the regime of magneti- 
zations and Lorentz factors we investigate. In Appendix 

fwe show how our results depend on the orientation of 
e field with respect to the simulation plane. 

3 In electron-ion shocks, the magnetization parameter as defined 
above pertains to ions. For electrons, <r e o = (mi/me) o" at injec- 
tion, where the electron kinetic energy is 7om e c 2 . So, it seems 
that the electron magnetization would depend on rrii/m e , for fixed 
a. However, on their way to the shock, electrons increase their 
average energy up to a fraction e e of the initial ion energy 70 m;c 2 , 
where e e is nearly insensitive to mi/m e . So, the "effective" electron 
magnetization <T e cf j- = (l/e e ) cr is insensitive to the mass ratio. 



3. ELECTRON-POSITRON SHOCKS 

In this section, we explore the physics of electron- 
positron shocks (rrii/m e = 1), focusing on the efficiency 
and rate of particle acceleration. In j ]3.1[ we discuss the 
long-term evolution of particl e ac celeration in unmagne- 
tized shocks (cr = 0), and in {3.2 we explore the physics 
of weakl y ma g netized flows, with 10~ 5 < a < 10" 1 . Fi- 
nally, in j ]3.3| we investigate the dependence of our results 
on the pre-shock bulk Lorentz factor 70 , for both unmag- 
netized and weakly magnetized shocks. 

3.1. Unmagnetized Shocks 



As shown by Spitkovsky ( 2008b ) and [ Sironi fc 
Spitkovskyl ( |2009b| ) in 2D and by [Spitkovskyf fl2005|) in 
3D, unmagnetized shocks are mediated by the so-called 



Weibel instability (|Weibel|195 9; Mcd vedev fc Loeb|1999 



Gruzinov fc Waxman|1999| |^| The free energy for the in 



stability comes from the counter-streaming between the 
incoming flow and a beam of high-energy particles re- 
flected back from the shock into the upstream ("return- 
ing" particles, from now on). The returning particles 
appear as a diffuse hot cloud to the right of the shock 
in the longitudinal phase space of Fig.Rjb), whereas the 
dense cold stream flowing into the shock is populated 
by the incoming plasma. The Weibel instability converts 
the kinetic energy of these two counter-propagating flows 
into small-scale (skin-depth) magnetic fields, organized 
in filaments stretching along the direction of propaga- 
tion of the shock (see the 2D plot of magnetic energy in 
Fig. ]4ta)). As the incoming flow approaches the shock, 
the filaments generated by the Weibel instability grow in 
strength and scale, eventually reaching sub-equipartition 
levels. At this point, they can efficiently randomize the 
bulk flow, and the shoc k forms. 
As demonstrated by Spitkovskyl (p008bl, unmagne- 



tized shocks in electron-positron plasmas are efficient 
particle accelerators, with ~ 1% of particles populat- 
ing a non-thermal power-law tail of slope p ~ 2. 40 In 
Fig. [2] we show the time evolution of the post-shock par- 
ticle spectrum, from the 2D simulation of an electron- 
positron unmagnetized shock. We confirm that the non- 
thermal population at late times can be described as a 
power-law tail of slope p ~ 2.4, and we show that the 
evolution of the high-energy part of the spectrum does 
not depend on the dimensionality of our computational 
domain. In fact, the inset of Fig. [2Fa) shows that the 
high-energy tail is the same in 2D (dotted line) and 3D 
(solid line). The difference at low energies, below the 
peak of the thermal Maxwellian that contains most of 
the particles, is just a consequence of the different adia- 
batic index between 2D and 3D, which changes the ana- 
lytic formula of a Maxwellian distribution (at low en- 
ergies, we have dNMB/dj oc 7 in 2D, as opposed to 

dNMB/d'y oc 7V7 2 — 1 in 3D). In summary, as regards 
to the efficiency of particle acceleration and the shape of 
the power-law tail, our results show that 2D simulations 
provide an excellent description of the 3D physics. 

4 More precisely, the instability results from the coupling of the 
electromagnetic filamentation (Weibel) instabil ity with the elec- 



trostatic two-stream instability, as explained by Bret (2009|. The 
resulting modes are oblique with respect to the streaming direction. 

5 We define the slope p of the non-thermal tail such that 
dN/d-y oc ~f~ p , where 7 is the particle Lorentz factor. 
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Fig. 2. — Left panel: Temporal evolution of the post-shock particle spectrum, from the 2D simulation of a 70 = 15 electron-positron 
unmagnetized shock. We follow the evolution of the shock from its birth (black curve) up to uipit = 15000 (red curve). The non-thermal 
tail approaches a power law with slope p = 2.4 (dashed line). In the subplot, we show at u) p \t = 1700 the comparison between 2D and 3D 
results (dotted and solid line, respectively). The difference at low energies is just a consequence of the different adiabatic index between 
2D and 3D. Right panel: The particle spectrum at different times is shifted along the rr-axis by (wpit/15000)- 1 / 2 and along the y-axis by 
(ojpit/lSOOO)^" 1 '/ 2 with p = 2.4, to show that the entire exponential cutoff scales in time as <x (wpit) 1 / 2 . The subplot shows that the 
maximum Lorentz factor 7max scales as oc (ojpit) 1 / 2 , both in 2D and in 3D (dotted and solid line, respectively). 



As shown in Fig. [2]ja), the high energy cutoff of the 
non-thermal tail of accelerated particles shifts in time to 
higher and higher energies (from black to red, as time 
progresses). The rate of evolution of the maximum par- 

'd), where 



re number 



ticle Lorentz factor j max is shown in Fig. [2'< 
Imax is defined as the Lorentz factor where t 
of particles drops by a factor of 10 5 with respect to the 
peak. Both in 2D (dotted) and in 3D (solid), the maxi- 
mum Lorentz factor grows in time as "/max oc (ojpii) 1 / 2 , 
much slower than the Bohm scaling (i.e., 7 max oc t) that 
is assumed in most models. 

The scaling ^ max oc t 1 ^ 2 is in agreement with analytical 
model s of particle diffusion in strong small-scale turbu 



lence ( |Shalchi k Dosch|2009||Plotnikoy et al |20n")|20i2| ), 

as appropriate for the fields generated by the Weibei in- 
stability in weakly magnetized shocks. In this regime, 
the isotropic diffusion coefficient for an accelerated par- 
ticle of energy e is 



D(s) ~ cA 



r L{e) 
A 



= cA 



eB\ 



(2) 



where A is the scale of the magnetic turbulence, B is the 
mean strength of the Weibel-generated fields, and r L (e) 
is the Larmor radius of a particle with energy e in the 



field B. As described by Achterberg et al. (|2001|) and 



Kirk k Reville (2010), the fact that the spatial diffu- 
sion coefficient scales as the square of the particle energy 
is simply a consequence of the small-scale nature of the 
Weibei fluctuations. The interaction of an accelerated 
particle with a Weibei filament of wavelength A results 
in a deflection of the particle momentum by an angle 
~ A/ri(e) <C 1. The expression in eq. ^ simply comes 
from the fact that it takes ~ [ri,(e)/X] 2 random scatter- 
ings in the upstream turbulence to diffuse the particle 
direction by roughly one radian, so that the particle can 
cross the shock back into the downstream and continue 



the acceleration process. Since the time needed to reach 
energy e in a relativistic shock is simply t ~ D /c 2 (apart 
from multiplicative factors of order unity) , the maximum 
particle Lorentz factor should scale in time as 



zB\ 



1/2 



(3) 



The strength B of the Weibel-generated fields can be 
parameterized in terms of the magnetic energy fraction 
e B = S 2 /87T7on i m.jC 2 , and the coherence length A of the 
Weibei fluctuations (measured in the direction transverse 
to the flow) can be written as a multiple A c /cj pi of the 
skin depth c/wpi, so that eq. (131 becomes 



7n 



7o 



(2esA c / (1Jpi a;pit) 



1/2 



(4) 



The scaling j max oc (wpii) 1 / 2 in our simulations is then 
expected if the combination egA c / u , is nearly constant 
in time. The temporal evolution of the magnetic struc- 
ture in unma gnetized shocks has been studied by |Keshet] 



et al. (20091). They found that the characteristic trans- 
verse scale of the upstream Weibei turbulence tends to in- 
crease with time, as particles are accelerated to higher en- 
ergies. However, the increase is not very significant, with 
A c / w growing by a factor of 2—3 between w p it ~ 10 3 and 

ujpit ~ 10 4 , corresponding to a tentative temporal scal- 
ing A c / Wp . oc (wpji) 1 / 4 , which is still sub-dominant with 
respect to the explicit temporal depende nce in eq. (El). 
As reg ards to the magnetic fraction €3, Keshet et al. 
(2009) found that 63 peaks at the shock, where its value 
is £B,sh — 0.1 at all times. We find that the magnetic en- 
ergy grows ahead of the shock with a longitudinal length 
scale increasing in time as Lb a f, while the particles 
are accelerated to higher energies. So, in a fixed volume 
around the shock, the average magnetic energy fraction 
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Fig. 3. — Temporal evolution of the Lorentz factor 7 of shock- 
accelerated particles (as a function of the time At acc spent in the 
acceleration process), from a sample of particles extracted from the 
2D simulation of an electron-positron unmagnctizcd shock. The 
particles are selected at Lo p \t = 4500 such that they belong to the 
non-thermal tail (i.e., 7 > 50 at ui pl t = 4500). Their mean Lorentz 
factor as a function of At acc is shown with the solid lines (red for 



positrons, blue for electrons), and it follows the <x At a 
indicated by the white dashed line. 



1/2 



scaling 



scales as eB.shLs ^ t- However, in the context of par- 
ticle acceleration, one should rather compute the mean 
magnetic energy within a volume defined by the diffusion 
length of the highest energy particles, that equals D/c 
and so scales as cx t. It follows that the mean magnetic 
energy fraction eg in a slab whose thickness is always 
the diffusion length of the highest energy particles is a 
constant fraction of 6b, sh, and so independent of time. 

In summary, the scaling ^max oc (ojpii) 1 / 2 is indeed 
expected in relativistic unmagnetized shocks mediated 
by the Weibel instability. From the inset of Fig. [^b) , we 
can measure the coefficient of the scaling, and we obtain 



It, 



To 



0.5(w pi i) 1/2 



(5) 



It is worth pointing out that the coefficient found above 
depends on our definition of "f ma x , taken here to be the 
Lorentz factor where the number of particles drops by a 
factor of 10 5 with respect to the peak (so, where jdN/dj 
is 10 5 times lower than the peak). For instance, we ob- 
tain a different coefficient (but the same temporal scal- 
ing) if jmax is defined by using 7 2 gL/V/g?7, rather than 
rydN/dj. The robustness of the scaling "/ ma x o^i^pit) 1 ^ 2 
is demonstrated in the right panel of Fig. ]2] where 
the particle spectrum at different times is shifted along 
the x-axis by (cjpit/15000)- 1 / 2 and along the y-axis by 
(wpii/lSOOO)^- 1 )/ 2 with p = 2.4. The fact that the expo- 
nential cutoffs of all the spectra overlap in the right panel 
of Fig. [^suggests that the downstream non-thermal tail 
can be aescribed at all times as a power law of fixed 
slope p = 2.4 terminating in an exponential cutoff that 
increases as cx (wpji) 1 / 2 . In other words, the scaling of 
the maximum energy as cx (wpji) 1 / 2 does not depend on 
our definition of "f m ax, provided that j ma x falls in the 
exponential cutoff of the particle distribution. 
The scaling with time in eq. ^ pertains to the highest 



energy particles accelerated in unmagnetized electron- 
positron shocks. In Fig. [3J we show that the same scal- 
ing describes the acceleration process of all the particles 
belonging to the non-thermal tail. We randomly select 
a sample of particles that at ui p it — 4500 have a Lorentz 
factor larger than 50 (i.e., they are part of the non- 
thermal component at that time), and we extract their 
energy histories at earlier times. For each particle, we 
compute the Lorentz factor 7 as a function of At acc , i.e, 
the time spent in the acceleration process (starting from 
the first crossing of the shock). Fig. [3] shows that, for 
both electrons (blue) and positrons (red), the mean en- 
ergy of the accelerated particles scales as cx (w p ; Atacc) 1 ^ 2 
(compare with the white dashed line, proportional to 
1/2 

Atacc)- In summary, by scattering off the small-scale 
Weibel turbulence, the energy of the particles injected 
into the acceleration process grows as cx {uJ v - l At acc ) 1 / 2 . 

In turn, this drives the evolution j ma x °c (wpji) 1 / 2 of the 
upper energy cutoff in the post-shock particle spectrum. 

3.2. Weakly Magnetized Shocks 

In this section, we investigate the physics and accel- 
eration properties of weakly m agnetiz ed shocks, in the 
regime 10 -5 < a < 10 _1 . In f 3.2.1 we show how the 
structure of the shock changes when varying the magne- 
tization, and in |3.2.2| we examine the dependence on a 
of the post-shock particle spectrum. 

3.2.1. Shock Structure 

Fig. [4] shows the flow structure from a set of 2D sim- 
ulations of electron-positron perpendicular shocks, with 
magnetization increasing from a = in the top row up 
to o~ = 10~ 2 in the bottom row. We plot the magnetic 
energy fraction €b — B 2 / Qir^QriimiC 2 in the left column 
(here B includes both the background ordered magnetic 
field B and the self-generated fields), and the positron 
longitudinal phase space x — ^iPi^x in the right column. 

With increasing magnetization from a = up to 
a = 10~ 2 (from top to bottom in Fig. [ij, the region 
around the shock filled with Weibel filaments becomes 
narrower, and for a > 10~ 2 the shock structure is en- 
tirely dominated by the shock-compression of the back- 
ground field Bq, with no evidence for Weibel-mediated 
turbulence. Here, the Larmor gyration of the incom- 
ing particles in the shock-compressed field is faster than 
the growth time of the filam entation instability, and th e 
Weibel modes cannot grow ( Lemoine &: Pelletier||2010[ ). 
This is in agreement with the results by SSOy and SSI1, 
showing that the Weibel instability is suppressed in per- 
pendicular relativistic shocks with a = 10 _1 . 

The right column of Fig. [4] shows that the longitudi- 
nal thickness of the upstream region filled with Weibel 
filaments is correlated with the typical distance traveled 
by the returning particles into the upstream (the return- 
ing particles populate the diffuse hot cloud ahead of the 
shock in the right column of Fig. [4J. The character- 
istic propagation length of the returning beam can be 
estimated as the Larmor radius of the returni ng parti 



cles in the backg round magnetic field B (see Lemoine 



& Pelletier 2010). As for the typical Lorentz factor ot 
the returning beam, we take the low-energy end of the 
non-thermal tail 7^ = JTinjTO; where rji n j ~ 5. Since 
the power-law slope of the non-thermal tail is p > 2, 
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Fig. 4. — Structure of the flow at u) p \t = 4500, from a set of 2D simulations of perpendicular electron-positron shocks with varying 
magnetization (from a = at the top to a = 10 -2 at the bottom, as indicated on the left of the figure). The left column shows the 2D plot 
of the magnetic energy fraction eg = B 2 / S-n^Qnirrnc 2 , and the right column shows the longitudinal phase space x — Tift x of positrons. 



the low-energy end will dominate both the number and 
the energy census. The longitudinal thickness of the up- 
stream region filled with Weibel turbulence will then be 



r Vinj I 



(6) 



It is worth pointing out that the magnetic energy profile 
will be desc ribed by this scale Lb, sat only at late times 
(see £ 3.2.2 for details on the particle spectrum at late 
times, and the role of Lg sat in setting the maximum 
energy of the shock-accelerated particles). For earlier 
times, the magnetic energy will grow in the upstream 
on a lon gitudinal scale evolving as Lb oc t, as described 
in [3.1 in the special case of a = shocks. Indeed, for 
cr = the expression in eq. ^ diverges, meaning that 
the characteristic width of the upstream region filled with 
Weibel filaments will steadily increase over time. 

As opposed to Lb, sat, which clearly depends on the 
magnetization of the flow (see the left column in Fig. El) , 
the magnetic energy at the shock always reaches ~ 10% 
of the kinetic energy of the pre-shock flow, i.e., tB,sh ~ 
0.1 regardless of a. Similarly, the transverse scale A of 
the Weibel filaments at the shock is nearly insensitive to 
a, and it equals A ~ 20c/a;pi (left column in Fig. H. 

As it is apparent in the magnetic energy plots oiFig. [4] 
(left column), the Weibel modes appear in 2D as mag- 



netic filaments stretched in the direction of propagation 
of the shock. Their 3D structure is shown in Fig. [5j 
for a relativistic electron-positron shock with magnetiza- 
tion <t = (top panel) and a = 10~ 3 (bottom panel). 
The background magnetic field Bq is oriented here along 
the z direction, in the same way as for our 2D simula- 
tions. The yz slice of the magnetic energy fraction in 
Fig. pic) shows that for a = 10~ 3 the magnetic field 
aheaa of the shock is primarily organized in pancakes 
stretched in the direction orthogonal to the background 
magnetic field (i.e., along y). This can be simply under- 
stood, considering that the Weibel instability is seeded 
by the focusing of counter-streaming particles into chan- 
nels of charge and current. In the absence of a back- 
ground magnetic field, the currents tend to be organized 
into cy lindrical filaments, as demonstrated by|Spitkovsky| 



1 2005 ) and shown in the yz slice of the top panel in Fig. 



In the presence of an ordered magnetic field along z, 
particles will preferentially move along the magnetic field 
(rather than orthogonal), so that their currents will more 
likely be focused at certain locations of constant z, into 
sheets elongated along the xy plane. This explains the 
structure of the magnetic turbulence ahead of the shock 
in the bottom panel of Fig. [5j common to all the cases 
of weakly magnetized shocks we have investigated (i.e., 
< a < 10" 1 ). 
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Fig. 5. — Structure of the flow, from the 3D simulation of an electron-positron shock with magnetization a = (top) or a = 10 -3 
(bottom). The xy slice shows the particle density (with color scale stretched for clarity), whereas the xz and yz slices show the magnetic 
energy fraction eg (with color scale stretched for clarity). 



3.2.2. Particle Spectrum and Acceleration 

We now explore the acceleration performance of weakly 
magnetized electron-positron shocks. In Fig. [6] we follow 
the evolution of the post-shock particle spectrum from 
the 3D simulation of a shock with magnetization a = 
1CT 3 . In Fig. m we compare the particle energy spectra 
at late times for different magnetizations, covering the 
range < a < 10 _1 . 

As compared to the results for unmagnetized shocks 
in Fig. [2j the evolution of the post-shock energy spec- 
trum in Fig. [6] for a flow with a — shows that the 
non-thermal tail initially grows to higher energies, but 
then it saturates (all the curves for io v \t > 1500 over- 
lap). The saturation of the high-energy tail is a robust 
result, holding in 3D (solid lines) and in 2D (dotted red 
line at uj p it — 3000), and it is clearly in contrast with 
the steady increase of the high-energy spectral cutoff ob- 
served for unmagnetized shocks in Fig. [2] The inset in 
Fig. [^confirms that the maximum Lorentz factor initially 
increases as j ma x oc (wpii) 1 / 2 (solid line for 3D, dotted 



for 2D), similarly to the case of unmagnetized shocks, 
but for ujpit > 1500 it saturates at j sat ~ 350. 

We find that the scaling ^ max cx (ojpit) 1 / 2 , followed by 
saturation at a constant 7 sa t, is a common by-product 
of the evolution of all magnetized relativistic shocks. In 
Fig. [7ja) we show several post-shock spectra for different 
magnetizations, after the non-thermal tail has reached 
the saturation stage. We cover the range 10 -4 < a < 
10 _1 , and for the sake of completeness we compare our 
results with the unmagnetized case a — 0, where the non- 
thermal tail is still evolving to higher and higher ener- 
gies. We find that strongly magnetized electron-positron 
shocks, with a > 10 -2 , are poor particle accelerators, 
in agreement with the conclusions of SS09. The post- 
shock spectrum at late times (see the black solid line for 
a = 10~ 2 ) is fully consistent with a Maxwellian distri- 
bution. This result does not depend on the reduced di- 
mensionality of our 2D computational domain. We have 
performed a large-scale 3D simulation of an electron- 
positron perpendicular shock with a — 10 , and we 
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Fig. 6. — Temporal evolution of the post-shock particle spectrum, 
from the 3D simulation of a 70 = 15 shock with <r = 10 -3 . We 
follow the evolution of the shock from its birth (black curve) up 
to UJpit = 3000 (red curve). The dotted red line is the particle 
spectrum at co p it = 3000 from a 2D simulation with the same pre- 
shock conditions. At late times, the particle spectrum reaches a 
steady state, and the maximum Lorentz factor saturates, as shown 
in the subplot (solid line for 3D, dotted for 2D). 



confirm that the post-shock particle spectrum (dotted 
cyan line in Fig. Fua)) does not show any evidence for 
particle acceleration. This undoubtedly proves that the 
absence of accelerated particles in the 2D simulations of 
perpendicular strongly magnetized shocks performed by 
SS09 is a physical consequence of the lack of sufficient 
self-generated turbulence j rather than an artifact of the 
re duced dimensionalit y ofthe simulation box, as argued 
by I Jones et~aT| ( |1998| ). 

For weakly magnetized shocks, with a < 3 x 10 -3 , we 
find efficient particle acceleration, with a non-thermal 



tail of slope p ~ 2.4 (dashed black line in Fig.Fna))) that 
contains ~ 1% of particles and ~ 10% of now energy, 
regardless of the magnetization. The low-energy end of 
the non-thermal tail does not significantly depend on the 
magnetization (7i„j ~ 570, or equivalently rji n j ~ 5), 
but the high-energy cutoff at saturation is systematically 
higher for lower magnetizations. This is confirmed by the 
inset of Fig. [7](a), where we plot the evolution in time of 
the maximum Lorentz factor ^max- Regardless of the 
magnetization, ~f max initially grows as j max oc (w pi i) 1/2 , 
with a coefficient of proportionality that does not signifi- 
cantly depend on a . At later times, the maximum energy 
departs from this scaling, and it saturates at a Lorentz 
factor 7 sa t which is larger for smaller magnetizations. 

For relatively high magnetizations (black for a = 10 -2 
and purple for a = 3 x 10~ 3 in the inset of Fig. (7^a)), 
the maximum energy initially grows as oc (wpit) 1 / 2 , then 
saturates at 7 sa i, and finally drops to a smaller value. 
The drop at late times is driven by a decrease in the 
strength of the Weibel turbule nce, as a result of h eating 
in the pre-shock medium (see Bret et al. 2010a for a 



6 More precisely, the fluctuations that get self-excited in <r > 
10 -2 shocks (cyclotron modes and their harmonics) have a short 
path length for emission and absorption, so they constantly enforce 
the local thermal equilibrium, giving Maxwellian energy spectra. 



discussion of temperature effects on the growth of beam- 
plasma instabilities). The pre-shock heating is induced 
by the so-called electromagnetic precursor wave, a train 
of transverse electromagnetic waves propagating into the 
upstream as a result o f the synchrotron maser instabil- 
ity at the shock front ( |Hoshino fc Arons 1991 Hoshino 
2008 1. The electromagnetic precursor heats the incom- 



ing particles, thus suppressing the Weibel instability and 
the particle acceleration at late times. Since the strength 
of the precursor wave de creases for smaller magnetiza- 



tions (Gallant et al. 



1992), for a < 10 3 the pre- heating 



induced by the electromagnetic precursor does not ap- 
preciably affect the structure of the magnetic turbulence 
and the process of particle acceleration. So, for a < 10~ 3 
the maximum Lorentz factor does not drop at late times. 

The scaling of the Lorentz factor at saturation j S at 
with respect to a is studied in the right panel of Fig. [7j 
In Fig. |7[d) we demonstrate that the magnetic energy 
profile ahead of the shock is characterized by a longi- 
tudinal scale Lb, sat oc a -1 / 2 , as suggested in eq. 
The curves in Fig. Qd) refer to oj p it = 6525, when 
shock has already reached a steady state, and the max- 
imum energy of the non-thermal tail has saturated (see 
Fig. [7]jc)). Since the magnetic turbulence that governs 
the acceleration process only extends across a region of 
longitudinal thickness Lb. sat ahead of the shock, the dif- 
fusion length of the highest energy particles will be lim- 
ited by the requirement D/c < is,satlj The Fermi pro- 
cess cannot proceed to higher energies, aue to the lack of 
sufficient magnetic turbulence on the relevant diffusive 
scales. More precisely, in Weibel-mediated shocks the 
non-thermal particles need to be continuously grazing 
the shock surface in order to be further accelerated, since 
the wavelength of the Weibel modes is too s mall to sus 
tain la rge-angle scatterings (as described by Spitkovsky 



2008b I . If a hypothetical particle were to penetrate into 
the upstream beyond the longitudinal scale L B , sa t of the 
Weibel turbulence, it would be deflected by the ordered 
background field away from its grazing trajectory into the 
downstream. In the absence of sufficient downstream tur- 
bulence, it will leave the shock transition region, prevent- 
ing further acceleration. The requirement D/c < Lb, sat 
sets the particle Lorentz factor at saturation to be 



Isat 
70 



1 /2 



(7) 



which implies, since e#, A c / Wp . and r\ in j do not signifi- 
cantly depend on a, that 7 sa t/7o oc a -1 / 4 . This scaling 
is demonstrated in Fig. Re)^ where we measure the co- 
efficient of proportionality tobe 



(8) 



Isat 4_ 

70 ~ C 1/4 

7 In reality, the Fermi-accelerated particles will sample both the 
upstream and the downstream turbulence, whereas Lg, sat only 
pertains to the upstream side of the shock. Since the decay length 
of the Weibel turbulence in the down stream is still a matter of 
debate l |Chang et al.|2008|[Keshet et al.|2009) , here we take Lg sat 
as our best estimate tor tne overall tmekness of the turbulent region 
around the shock, including both upstream and downstream. 

8 In Fig. [7je), the open circles for a = 3 X 10~ 3 and a = 10 -2 
show the Lorentz factor that is reached during the plateau phase 
after the initial growth as oc (w^t) 1 / 2 , and before the drop at late 
times induced by the electromagnetic procursor (see the black and 
purple lines in Fig. |7jc)). 
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Fig. 7. — Left panel: Dependence of the post-shock particle spectrum on the upstream magnetization, from a set of 2D simulations of 
electron-positron shocks with 70 = 15. We vary the magnetization between a = (red curve) up to a = 10 — 2 (black curve), showing that 
the Fermi process is suppressed with strong pre-shock fields. This is confirmed by the post-shock spectrum of a 3D simulation with a = 10 _1 
(dotted cyan line). In the inset, we follow the maximum particle Lorentz factor over time, from the growth as "/max oc (uipit) 1 ^ 2 to the 
saturation at a constant "f sa t- Right panel: The particle spectrum for different magnetizations is shifted along the x-axis by (a-/10 -4 ) 1//4 
and along the j/-axis by (10~ 4 /<r)( p ~ with p = 2.4, to show that the entire exponential cutoff scales with magnetization as oc tr -1 / 4 . 
In the inset (d), we show that the magnetic energy profile ahead of the shock has a characteristic longitudinal scale Lb, sat oc cr -1 / 2 . In 



the inset (e), we show that the Lorentz factor at saturation -jsat scales with the magnetization as -y S at oc a 



-1/4 



Finally, the right panel of Fig. [7] (yellow curve for a = 
10~ 4 , green for 3 x 10~ 4 and blue for 10~ 3 ), where 
we rescale the x-axis by (er/lO -4 ) 1 / 4 and the y-axis by 
(lCr 4 /cr) (p ~ 1)/4 with p = 2.4, demonstrates that the en- 
tire exponential cutoff of the non-thermal tail scales as 
oc ct -1 / 4 , and so the proportionality j S at/lo oc er~ 1//4 
does not depend on the details of our definition of 7 sa t 
(yet, the coefficient of the scaling in eq. ^ does depend 
on the definition of J sa t)- 

3.3. Dependence on the Lorentz Factor 70 

In Figs. [8] and [9j we show the dependence of the post- 
shock particle spectrum on the upstream bulk Lorentz 
factor 70, for both unmagnetized (Fig. [8| and weakly 
magnetized (er = 10 -3 ; Fig. [9| electron-positron shocks. 
We vary 70 from 70 — 3 (a mildly rclativistic shock) up 
to 70 = 150 (an ultra-relativistic shock). 

For both unmagnetized and weakly magnetized shocks, 
we find that the acceleration physics does not depend on 
the bulk Lorentz factor of the upstream flow, if 70 > 10. 
Both the acceleration efficiency and the rate of particle 
energization are insensitive to 70, in the limit of ultra- 
relativistic flows. In particular, the insets in Figs. [8] and 
[9] show that the temporal evolution of "/max/ 7o is the 
same for all 70 > 10, as predicted by eqs. ([5]) and pi). As 
shown in the inset of Fig. [9J for a — 10~ 3 both the scaling 
7max oc (wpit) 1 / 2 at early times and the saturation of the 
maximum energy at late times hold regardless of 70 > 10. 

The decrease in the acceleration efficiency for 70 < 10 
is a consequence of the fact that the Weibel and oblique 
instabilities are suppressed when the beam of returning 
particles has a significant transverse dispersion (i.e., it 
is not cold, in the pre-shock frame). The instability will 
be quenched if u g < k ■ Ad, where oj g is the growth 
rate of the most unstable mode, that usually occurs for 



k ~ k± ~ Wpi/cQ In the post-shock frame, the parallel 
and perpendicular components of the momentum of the 
returning beam are pu ~ p± ~ r yinj m i c - I n t ne pre-shock 
frame, pj, ~ 7of|| and p'^ ~ p±, so that the charac- 
teristic beam transverse velocity in the pre-shock frame 
is Av'j_ ~ cpj_/p(i ~ c/70, which is larger for smaller 70. 
The beam transverse dispersion should be com pared with 



the gro wth rate of the relevant instability (e.g.,|Bret et al. 



2010b I, which in the upstream frame of electron-positron 
shocks is oj 



for the Weibel mode and 



~ (£& / '77 •mj) 1 ' ' 3 'w pi for the oblique mode. Here, £b is the 
density ratio between the beam of returning particles and 
the incoming flow, measured in the downstream frame, 
whereas r)i n j = 7mj/7o- Both and r\i n j are nearly in- 
sensitive to 70 (Figs. [8] and [9] show that the non-thermal 
tails for different 70 are all starting at the same Ji n j/jo> 
and with similar normalizations), so that the generation 
of Weibel and oblique modes in electron-positron flows 
should be suppressed for 



7o < (6/W 1/2 [Weibel] 
7o 



< (6/%ni)" V3 [Oblique] 



(9) 
(10) 



In turn, the suppression of the Weibel and oblique in- 
stabilities, which are responsible for the generation of 
the magnetic turbulence that governs the Fermi process, 
results in the poorer acceleration capabilities of mildly 
relativistic shocks, as compared to their ultra-relativistic 
counterparts (see the green and yellow curves in Figs. [8] 
and|9j for 70 = 5 and 70 = 3 respectively). 



9 Here we only consider electron-positron plasmas, so c/uj p \ 
c/u) pc . For electron-ion flows, see |Lemoine fc Pelletier| | |2011^ . 
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Fig. 8. — Dependence of the post-shock particle spectrum on the 
upstream bulk Lorentz factor 70, from a set of 2D simulations of 
unmagnetized electron-positron shocks. The spectra are shifted 
along the x-axis by 15/70, to facilitate comparison with the ref- 
erence case 70 = 15. In the subplot, we show that the maximum 
Lorentz factor 7 ma x scales as <x (wpit) 1 / 2 for all the values of 70 we 
explore. Both the main plot and the subplot show that the physics 
of relativistic shocks does not depend on 70, for 70 > 10. 
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Fig. 9. — Dependence of the post-shock particle spectrum on the 
upstream bulk Lorentz factor 70, from a set of 2D simulations of 
cr = 10 -3 electron-positron shocks. See the caption of Fig. [8] for 
further details. For all the values of 70, the maximum particle 
energy saturates after u) p it ~ 1500. 



4. ELECTRON-ION SHOCKS 

In this section, we investigate the physics of electron- 
ion shocks, by employing a reduced mass ratio Wi/m e = 
25 or rrii/m e — 100. As we describe in Appendix |A[ our 
conclusions do not change for higher mass ratios, up to 
the realistic value irii/m e ~ 1836. In |4.1| we show how 
the structure of the flow changes for diffe rent magneti- 
zations, in the regime < a < 10~ 3 . In \ A.2 we follow 
the time evolution of the post-shock spectrum, and we 
comment on the effect of the flow magnetization on the 
long term evolution of the acceleration physics. Finally, 



bulk Lorentz factor of the pre-shock flow. 

4.1. Shock Structure 

The physics of elec tron-ion unmagnetiz ed sh ocks has 
been investigated by [Spitkovsky ( 2008a) and Martins| 
et al.| ( |2009[ ) in 2D and by |Maugb0lle| ( ]2011h in 3D. They 
have found that electron-ion unmagnetized (i.e., a = 0) 
shocks are mediated by the Weibel instability, in the 
same way as electron-positron shocks. The instability is 
seeded by the counter-streaming between the incoming 
flow and the ions and electrons reflected back from the 
shock into the upstream (see the diffuse hot cloud of ions 
and electrons moving ahead of the shock in Fig.fTolb) and 
(c), respectively). As they approach the shocK^the in- 
comin g electrons are heated by the time- varying W eibel 



fields (Gedalin et al. 2012 Plotnikov et al. 2012 1, and 



they reach equipartixion with the 10ns before entering 
the shock. As the electrons get heated on their way to 
the shock, the typical transverse scale of the Weibel fil- 
aments increases from the electron skin depth c/w pe far 
ahead of the shock to the proton skin depth c/w p i just in 
front of the shock. In the 2D plot of the magnetic energy 
in Fig. 10 ^a), the transverse scale of the Weibel filaments 
at the sh ock amounts to a few ion skin depths. 

In Fig. [lOj we show how the 2D structure of the mag- 
netic energy (first column) and the longitudinal phase 
spaces of ions and electrons (second and third columns, 
respectively) change for increasing magnetization, as 
compared to the unmagnetized case (top row) most com- 
monly discussed in the literature. We investi gate the 



10- 



in f 4.3 we explore the dependence of our results on the 



range < a < 10~ 3 , from top to bottom in Fig. 

Similarly to the case a — 0, in the regime a 
of weakly magnetized shocks the flow structure is dom- 
inated by the Weibel modes. In the same way as for 
electron-positron magnetized flows in §3.2. 1[ the thick- 
ness of the upstream region that is filled with Weibel 
filaments is sensitive to the flow magnetization, being 
set by the typical Larmor radius of returning ions, that 
scales as oc cr" 1 / 2 (see eq. The two main differ- 

ences with respect to electron-positron shocks are in the 
transverse scale of the Weibel filaments and in their ori- 
entation relative to the shock normal, as we now discuss. 

As the magnetization increases in the range 10~ 5 < 
a < 10 -4 , the Weibel filaments ahead of the shock ap- 
pear narrower. As anticipated above, if the characteristic 
electron thermal Lorentz factor in the pre-shock frame 
is 7tfc,e; the fastest growing Weibel modes will operate 
on a transverse scale ~ ^/lth,e c/oj pc . For a < 10~ 5 , 
the region filled with Weibel filaments is wide enough 
to drive efficient ion-to-electron energy transfer ahead of 
the shock, and jth,e ~ rni/m e as the electrons enter the 
shock. It follows that the characteristic transverse scale 
of the Weibel filaments in Fig. [To| d) is comparable to the 
ion skin depth c/w p i = \J mi/m e c/ui pc for a < 10 -5 . On 
the other hand, for 3 x 10~ 5 < a < 10~ 4 , the electrons 
ahead of the shock do not have enough time to reach 
equipartition with the ions, in the Weibel turbulence gen- 
erated ahead of the shock. So, 7th, e stays smaller than 
rrii/m e , and the Weibel filaments are narrower. 

Since for 3 x 10~ 5 < a < 10~ 4 the electrons enter 
the shock with energy lower than the ions (resulting in a 
lower dispersion of their post-shock longitudinal phase 
space, compare panels (k) and (1) for a = 10 -4 ), a 
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Fig. 10. — Structure of the flow at LO^t = 2250, from a set of 2D simulations of perpendicular electron-ion shocks (m,i/m e = 25) with 
varying magnetization (from a = at the top to <r = 10~ 3 at the bottom, as indicated on the left of the figure). The first column shows 
the 2D plot of the magnetic energy fraction eg = B 2 / '87T7on;miC 2 , the second column the longitudinal phase space x — Tjft x of ions, and 
the third column the longitudinal phase space x — "f e /3 e ,x of electrons. 



smaller fraction of electrons will be injected into the ac- 
celeration process, relative to ions. This is demonstrated 
in the longitudinal phase spaces of a = 3 x 10~ 5 (panel 
(h) for ions and (i) for electrons) and a = 10~ 4 (panels 
(k) and (1)), that show a smaller number of returning 
electrons, with respect to ions. By comparison, in the 
cases where the ion-to-electron energy transfer ahead of 
the shock is efficient (i.e., a < 10~ 5 ), returning elec- 
trons and ions have comparable number density and en- 
ergy. If ions outnumber electrons, as it is the case for 
3 x 1CT 5 < a < 1CT 4 , the orientation of the Weibel 
modes ahead of the shock will be determined by the di- 
rection of the electric current seeded by the returning 
ions. For a background magnetic field in the +z direc- 
tion, as we employ here, the ions back-streaming ahead 
of the shock will preferentially move with p y ~ —p x < 0. 
The electric current they generate is then in the same 
direction as the Wcibcl filaments seen in the upstream 
region of a = 3 x 10~ 5 (panel (g)) and a = 10~ 4 (panel 
(j)) shocks. In contrast, if the injection of ions and elec- 
trons into the acceleration process is equally efficient (so 
that the two species equally contribute to the popula- 
tion of returning particles), the mean flux of returning 
particles stays aligned with the shock normal along +x, 
and the Weibel filaments are stretched along the shock 
direction of propagation, as shown in panel (a) for a = 
and panel (d) for a = 10~ 5 . 

For highe r m agnetizations (a = 10~ 3 in the bottom 
row of Fig. 10), the returning ions are confined closer 
to the shockoy the background magnetic field (panel 
(n)), no returning electrons are present (panel (o)), and 
the structure of the shock is no longer dominated by 
the Weibel filaments (panel (m)). The 2D plot of mag- 



netic energy in Fig. 10 ^m) shows that the shock surface 
is rippled, on a scale comparable to the ion Larmor ra- 
dius. The phenomenon of shock ripplin g is well known 
in non-relativistic supercritical shocks (Burgess 20061, 
where it is associated with the presence of gyrating re- 
flected ions in the shock transition layer. In our simula- 
tions of relativistic perpendicular shocks, we observe the 
shock ripples only in a limited regime of magnetizations, 
3 x 10 < a < 10 -1 . For both smaller and higher mag- 
netizations, the rippling is suppressed by the pre-shock 
heating of the incoming electrons. The cause for the 
electron heating varies with the flow magnetization (the 
time-dependent Wcibcl fields for 3 x a < 10~ 4 ; the elec- 
tromagnetic precursor wave for a > 10 , as discussed 
by SS11), but in all the cases it results i n a suppression 
of the shock ripples. As we discuss in SfOl we do not 



find that shock rippling facilitates the injection of parti- 
cles into the acceleration process. This result has been 
confirm ed with a simulation box three times as large as 
in Fig. fl0[m), such that to accommodate three shock 



ripples m the computational domain. 

4.2. Particle Spectrum and Acceleration 

We now investigate the acceleration physics in perpen- 
dicular electron-ion magnetized shocks, by employing a 
reduced mass ratio mi/m e — 25. As we show in Ap- 
pendix [Aj our conclusions do not change for higher mass 
ratios, u p to the realistic value nii/m e ~ 1836. 

In Fig. [IT] we follow the evolution of the post-shock par- 
ticle spectrum from the 2D simulation of a perpendicular 
shock propagating in a weakly magnetized plasma, with 
magnetization a = 10 -5 . Similarly to the case of unmag- 
netized electron-ion shocks (i.e., a = 0) investigated by 
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Fig. 11. — Temporal evolution of the post-shock particle spec- 
trum, from the 2D simulation of a 70 = 15 electron-ion (m,i/m e = 
25) shock in a flow with magnetization a = 10 — 5 . We follow the 
evolution of the shock from its birth (black curve) up to ui p it = 2500 
(red curve). In the top panel we show the ion spectrum, in the 
bottom panel the electron spectrum. The ion non-thermal tail ap- 
proaches at late times a power law with slope p = 3.0, whereas 
p = 2.5 for electrons (dashed lines). In the bottom panel, we over- 
plot the ion spectrum at co p it = 2500 with a red dotted line. Inset 
(a): mean post-shock ion (red) and electron (blue) energy, in units 
of the bulk energy of the upstream flow. The dashed blue line 
shows the electron energy at injection. Inset (b): temporal evo- 
lution of the maximum Lorentz factor of ions (red) and electrons 
(blue), scaling as <x (cjpii) 1 / 2 at late times (black dashed line). 



ISpitkovsky l ([ 2008a) an d [Martins et al.| ( |2009|) in 2D and 
by |Hau gb0llc ( 2011| ) in 3D, both the ion and the electron 
spectra show a prominent non-thermal power-law tail of 
accelerated particles, beyond the thermal distribution. 

As a result of efficient energy transfer in front of the 
shock, mediated by the Weibel turbulence, the electrons 
enter the shock nearly in cquipartition with the ions. 
This is demonstrated in the inset (a) of Fig. 11 where we 
follow the temporal evolution of the mean electron and 
ion energies in the post-shock flow, in units of the initial 
ion energy. At late times, the mean electron energy (blue 
line) is a fraction ~ 40% of the initial ion bulk energy, and 
it amounts to a fraction ~ 65% of the mean post-shock 
ion energy (red line). This is much larger than the initial 
electron energy in the pre-shock flow, which is shown as 
a dashed blue line in the ins et (a) for comparison. In 
the bottom panel of Fig. 11 the energy equipartition 
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Fig. 12. — Dependence of the post-shock particle spectrum on the 
upstream magnetization, from a set of 2D simulations of electron- 
ion shocks (mi/m e = 25) with 70 = 15. In the top panel we show 
the ion spectrum, in the bottom panel the electron spectrum. We 
vary the magnetization between a = (red curve) up to a = 10 — 3 
(black curve), showing that the Fermi process is suppressed for 
large a. In the insets, we follow the maximum particle Lorentz 
factor over time, for ions (panel (a)) and electrons (panel (b)). 



between ions and electrons at uj pi t = 2500 explains the 
proximity of the thermal peaks in the electron (red solid 
line) and ion (red dotted line) spectra. 

Since electrons and ions enter the shock with nearly 
the same energy, the injection efficiency into the acceler- 
ation process will be similar for the two species. In fact, 
by comparing the red solid line (for electrons) and the red 
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dotted line (for ions) in the bottom panel of Fig. 
show that the non-thermal tails in the electron and ion 
spectra have comparable normalizations, yielding an ac- 
celeration efficiency for both species of roughly ~ 1% by 
number and ~ 10% by energy. The power law in the ion 
non-thermal tail is steeper than for electrons (p ~ 3.0 for 
ions and p ~ 2.5 for electrons), but it may evolve toward 
flatter slopes at later times, so that for both species the 
spectral index will asymptote to the value p ~ 2.5 found 
in electron-positron shocks. 

For both ions and electrons, the non-thermal tail ex- 
tends in time to higher and higher energies. As shown 
in the inset (b) of Fig. 11 the maximum energy at late 

x / 2 , for both ions (red line) and 



times increases as oc (w p it) 
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electrons (blue line) This is the same scaling observed 
in Fig. [2] for electron-positron shocks. The similarity be- 
tween electron-positron and electron-ion shocks is to be 
expected, since in electron-ion flows the two species enter 
the shock nearly in energy equipartition, so our electron- 
ion shocks should display the same acceleration physics 
present in electron-positron flows. From the inset (b), 
we find that at late times (u> p it > 700) 
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0.25 (w pi i) 1/2 



(11) 



Even though the scaling of the maximum energy as oc 
7o(w B j£) 1 / 2 is the same as in electron-positron shocks (see 
eq. (pi)), the coefficient of proportionality in electron-ion 
shocks is smaller by a factor of ~ 2. This simply reflects 
the position of the thermal peak, located at j p k ~ 70 
in electron-positron flows and at i P k,i ~ 7p*,eW& e /mj ~ 
7o/2 in electron-ion flows, just as a result of efficient ion- 
to-electron energy transfer ahead of the shock (see the 
location of the thermal peaks in Fig. [IT] as compared to 
the bulk Lorentz factor 70 = 15). 

With a similar argument, we expect that in electron- 
ion flows with magnetization a the maximum Lorentz 
factor should saturate at 
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(12) 



which replaces eq. ([8]), valid for electron-positron plas- 
mas. For a — 10 , the expected saturation value would 
be Jsat,i ~ lsat,e m e/'m-i — 550, well beyond the max- 
imum Lorentz factor reached at the end of our simula- 
tion. We would need to evolve the simulation up to much 
longer times (w p it ~ 18000) to capture the saturation of 
the Lorentz factor in our a = 10 -5 electron-ion shock. 

We then extend our investigation to differ ent magneti- 
zations, in the range < a < 10~ 3 . In Fig. n2l we show 
how the post-shock spectrum of ions (top panel) and elec- 
trons (bottom panel) depend on the magnetization of the 
flow. We find that shocks with a < 3 x 10~ 5 are efficient 
particle accelerators, and the spectrum shows a promi- 
nent non-thermal tail containing ~ 1% of particles and 
~ 10% of flow energy (for both ions and electrons). The 
exponential cutoff of the non-thermal tail extends in time 
to higher energies, and at late times it follows the scaling 
oc (ajpii) 1 / 2 discussed above, for both ions (inset (a)) and 
electrons (inset (b)). The coefficient of proportionality 
tends to be lower for higher a (compare the green line, 
for a — 3 x 10~ 5 , to the red line, for a = 0). This minor 
dependence on a might be caused by the slight decrease 
in the transverse wavelength A of the Weibel modes with 
increasing magnetization, as apparent in Fig. 10 (com- 
pare panels (a), (d) and (g) for a = 0, a = 10~ 5 and 
a = 3 x 10~ 5 , respectively). In turn, the transverse scale 
A of the Weibel fluctuations enters the temporal scaling 
of the maximum energy, as illustrated in eq. ([4]). How- 
ever, the effect of the magnetization on the scaling of the 
maximum energy is of minor importance, and we will 
neglect it in the following. 



10 |Martins et ah] l |2009| reported a different scaling for the tem- 
poral evolution ot the energy of shock-accelerated particles, but 
their conclusions might have been misguided by the limited times- 
pan of their simulation runs. 



With increasing magnetization, the acceleration effi- 
ciency decreases (blue lines for a — 10~ 4 ), and for 
a > 10~ 3 the post-shock spectrum does not show any 
evidence for non-thermal particles. For a — 10~ 3 (black 
lines) , the downstream electrons populate a thermal dis- 
tribution, that peaks at a lower energy than for smaller 
magnetizations, since the Weibel modes that mediate 
the electron heating ahead of a < 10~ 4 shocks are now 
suppressed. The non-Maxwcllian shape of the ion spec- 
trum for a = 10~ 3 results from incomplete ion thermal- 
ization at the shock front, and it will eventually relax 
to a Maxwellian distribution further downstream from 
the shock. The absence of non-thermal particles for 
a = 10~ 3 proves that shock rippling in electron-ion rela- 
tivistic flows (see Fig. |10[ m)) is not a promising mecha- 
nism to mediate the injection of particl es in to the Fermi 
process. In the insets (a) and (b) of Fig. 12 the suppres- 
sion of non-thermal particle acceleration in a > 10~ 4 
shocks is revealed by the fact that the maximum Lorentz 
factor of ions (panel (a)) and electrons (panel (b)) is 
nearly constant with time (blue and black lines). 

Finally, we refer to SS11 for a detailed investigation of 
the acceleration physics of electron-ion shocks at higher 
magnetizations (a > 10~ 2 ), a regime where the electro- 
magnetic precurs or wave resulting from the synchrotron 
maser instability ( jHoshino fc Arons|199l| |Hoshino|2008| ) 
can appreciably affect the pre-shock flow. 

4.3. Dependence on the Lorentz Factor 70 

Fig. [13] shows how the spectrum of electron-ion 
(mi/m e = 100) magnetized shocks depends on the bulk 
Lorentz factor 70, by comparing the cases 70 = 60 and 
240 with our reference value 70 = 15. For both ions 
(top panel) and electrons (bottom panel), we find that 
our results are basically the same regardless of_7o , mod- 
ulo an overall shift in the energy scale. Fig. |13| confirms 
that shocks with a = 10~ 4 are poor particle accelerators 
(blue lines), whereas the generation of a prominent non- 
thermal tail is observed for a < 3 x 10 -5 (green line for 
a = 3 x 10~ 5 , yellow line for a = 10~ 5 ). 

Most importantly, the spectra in Fig. [13] obtained for a 
mass ratio mj/m e = 100, are entirely consistent with the 
results in Fig. [l2| referring to mi/m e = 25. This suggests 
that the acceleration physics in electron-ion shocks can 
be confidently captured with a reduced mass ratio, as we 
will further demonstrate in Appendix \K\ 

5. ASTROPHYSICAL APPLICATIONS 

In this section, we explore the implications of our re- 
sults for the accelera tion process at the external shocks of 
GRB afterglows (in f 5.1 1 and at t he te rmination shock of 
Pulsar Wind Nebulae [FWNe; in f 5.2 ). Our discussion of 
the termination shock in pulsar winds might equally ap- 
ply to the termination hot spots of AGN jets - f or exam- 
ple, the par ameters of the Cygnus A hot spots ( |Stawarzj 
et al. 2007[ ) are surprisingly similar to the termination 
shock of the Crab Nebula. 

We make use of the results presented in the previous 
sections. In particular, we employ the fact that the evo- 
lution in time of the Lorentz factor of shock-accelerated 
particles follows the scalings 
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Fig. 13. — Dependence of the post-shock particle spectrum on the 
upstream bulk Lorentz factor 70, from a set of 2D simulations of 
electron-ion (m;/m e = 100) shocks with different magnetizations 
(yellow for a = 10~ 5 , green for a = 3 X 10 — 5 , blue for a = 10 -4 ). 
In the top panel we show the ion spectrum, in the bottom panel 
the electron spectrum. 
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0.25 7o ( Wpi i) 1 / 2 (14) 



in electron-positron and electron-ion shocks, respectively. 
Here, 70 is the bulk Lorentz factor of the upstream flow 
(in the post-shock frame) and w p i indicates either the 
positron plasma frequency (for electron-positron flows) 
or the ion plasma frequency (for electron- ion shocks). 

The increase of the maximum particle energy over time 
proceeds up to a saturation Lorentz factor that is con- 
strained by the magnetization a of the upstream flow 
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(15) 
(16) 



in electron-positron and electron-ion shocks, respectively. 
Further energization is prevented by the fact that the 
self-generated turbulence is confined within a region of 
thickness Lb, sat ~ ct -1 / 2 c/w P i around the shock. 

5.1. GRB Afterglows 

The afterglow emission in GRBs is usually attributed 
to synchrotron radiation from electrons accelerated at 
the forward shock (the so-called "jitter" paradigm does 
not seem to relevant in GRB afterglows, as shown by 



Sironi & Spitkovsky 2009b ) . The shock propagates either 
in the ISM, with constant density, or in the stellar wind of 
the GRB progenitor, with density scaling in radius as n oc 
R~ 2 . As discussed in sjl] the magnetization of the ISM 
is a ~ 10~ 9 , and it is independent of the shock radius 
R. The same holds for a stellar wind, if the magnetic 
field in the wind is primarily toroidal, thus decreasing 
as oc R^ 1 . We ass ume that, as appropriate for a Wolf- 
Raye t progenitor (Eichler & Usov 1993 Sagi & Nakar 



20121, the number density in the wind at R = 10 iB cm 



is n = 0.3 Ti-o.5 cm j an d the strength of the magnetic 
field at R = 10 18 cm is B w = 10~ 5 B W _ 5 G, so that the 
magnetization parameter in the wind will be 



4mnmjC 2 



1.7 x \Q- % B 2 



w. 



-1 
-0.5 



(17) 



We treat separ ately t he case of a blast wave pr opaga ting 
into the ISM (j j5.1.1| ) and into a stellar wind (j ]5.1.2| . 

In the following, we study the process of particle accel- 
eration in GRB afterglows by assuming that the shock 
properties are only determined by the instantaneous 
magnetization and Lorentz factor of the flow. In other 
words, we neglect the feedback of particles and magnetic 
turbulence produced at early times on the shock struc- 
ture and the acceleration process at later times. 

Also, we only focus on the particles reaching the high- 
est energies, with the goal of assessing whether GRB af- 
terglow shocks are promising candidates for the accel- 
eration of the UHE CR protons obse rved by the Pierre 
Auger Observatory (Abraham 2010[ ), and whether syn- 
chrotron radiation from the shock-accelerated electrons 
can explain the early sub-GeV emission of Fermi-LAT 
GRBs (e . 
lellini et al 



De Pasquale 2010 Ghis 



Ac kermann 2010 

lUb- Our st rategy parallels closely tiie study 
by I Sagi fc JNakar] ( |2012[ ), with one important difference. 
They were inferring the rate of electron acceleration in 
GRB external shocks from the detection of high-energy 
photons at early times, whereas our PIC simulations can 
provide from first principles an estimate for the accelera- 
tion rate. Our goal is then to verify whether synchrotron 
emission from the accelerated electrons can explain the 
highest energy photons detected by Fermi-LAT. 

Since we only focus on the highest energy electrons, 
we do not attempt to provide a complete description of 
the full spectral and temporal evolution of GRB after- 
glows. For this, we would need a self-consistent estimate 
for the overall extent of the region populated with Weibel 
turbulence. A strict lower limit is provided by the thick- 
ness Lb, sat of the upstream layer filled with Weibel fila- 
ments (see eq. ([6])), but the decay length of the Weibel- 
generated fiel ds m the downstream region is still a m at- 
ter of debate ( |Chang et aL|2008| |Keshet et aL|2009l) . A 
physically-grounded estimate for the overall extent of the 
region with Weibel turbulence could be used to predict 
the temp oral and spectral evo lutio n of GRB afterg lows, 
following |Rossi fc Rees| ( |2003[ ) and |Lemoine| ( |2012| ) . 



5.1.1. ISM 

Under the assumption of adiabatic expansion in a con- 
stant density medium, the bulk Lorentz factor of the 
shock during the relativistic deceleration phase will be 
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described by the |Blandford & McKee| (|1976| solution 

(18) 
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at distances larger than the deceleration radius Rdec = 
(nEo/WirrlnmiC 2 ) 1 / 3 . Here, E is the isotropic- 
equivalent explosion energy, Tq is the initial Lorentz fac- 
tor of the blast wave, and n is the ISM number density 
(measured in the rest frame of the ISM). By defining 



E = 10 54 £; 0; 54 ergs, T = 10 2 - 5 r 0i2 . 5 and n 
the deceleration radius amounts to 



no cm 



Rde 



1.6 X 1U £'o,54 i 0,2.5 n Cm 



(19) 
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The shock becomes non-relativistic 
{nE a /l&i:nm l c 2 ) 1 /' i 

which point our results will no longer be applicable 

We now evaluate the acceleration performance of GRB 
external shocks, first for protons and then for electrons. 
We use the fact that the Lorentz factor 70 we have em- 
ployed in the previous sections coincides with the in- 
stantaneous Lorentz factor of the GRB blast wave, i.e., 
70 = r at all radii. The maximum Lorentz factor of the 
shock-accelerated protons is constrained by the magne- 
tization of the pre-shock flow, as illustrated in eq. (16). 
When transforming into the upstream fr ame (coincindent 
with the ISM frame), the result in eq. (16) needs to be 



multiplied by an additional factor of T, giving 



up 
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8.0 x 10 7 #o,54 1 V-l /iR i7 > (2°) 



where a = 10 -9 <7_g and R = 10 17 i?i7 cm. This strict 
upper limit could be circumvented if the coherence scale 
Aco/i of the ISM magnetic field in the longitudinal direc- 
tion (i.e., along the shock normal) is small enough such 
that the acceleration of protons up to 7"^ i takes longer 
than the passage of the shock through a region of length 
X co h- From the acceleration timescale in eq. ( |14[ ), this ar- 
gument gives an upper limit on the required longitudinal 
coherence length of the ISM field (in the ISM frame) 

X coh < 2.2 x 10 16 £ V 5 W ^-l /2 Ri7 /2 cm , (21) 



where in eq. ( 14 ) we have used the fact that the proton 



plasma frequency in the ISM is 



-'pi 



1.3 x 10 3 rin /2 Hz 



(22) 

The coherence scale of the ISM magnetic field is believed 



to be a few parsecs, so the constraint in eq. (21) is not 
likely to be satisfied, and proton acceleration to Lorentz 
factors larger than 7"^ i will not occur in the ISM. 

A separate constraint comes from the re quir ement that 
the proton acceleration timescale in eq. ( 14 ) should be 
shorter than the age of the blast wave ~ 77/Tc as mea- 
sured in the post-shock frame. This constrains the max- 
imum proton Lorentz factor to be, in the ISM frame, 
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(23) 



By comparing eq. ( 20 1 and eq. ( 23 ), we see that the max 



imum proton Lorentz factor 7^^ = min[^ t>i , lT ge ,i\ 
is always limited by the mag netization of the pre-shock 
flow, as prescribed in eq. (20). From this, it is clear that 



the maximum energy of protons accelerated in GRB ex- 
ternal shocks propagating into the ISM is too small to 
explain the extreme Lore ntz factors (~ 10 11 ) of the UHE- 
CRs observed by Auger ( Abraham||2010 ). 

We now turn to the acceleration ot electrons. The 
Lorentz factor of the thermal peak in the electron dis- 
tribution will be, in the post-shock frame, 
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(24) 



where the factor of two in the denominator comes from 
the assumption of energy equipartition between electrons 
and protons behind the shock. The Lorentz factor "f p k, e 
is a strict lower limit for the low-energy end of the non- 
thermal tail (or equivalcntly, for the minimum Lorentz 
factor 7 rra i I i j e of the shock-accelerated electrons), so we 
suggest to take 7 mi „, e ~ 7 p a. 



With regards to the acceleration of electrons, eq. ( 16 1 
gives an absolute upper limit constrained by the now 
magnetization, which in the post-shock frame yields 
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(25) 



As explained above, this strict upper limit implicitly as- 
sumes that the ISM field does not change its orienta- 
tion as electrons are accelerated up to the Lorentz factor 
7sat,e, which is generally expected given the parsec-scale 
coherence length of the ISM magnetic field. 

Radiative energy losses are likely to limit the maximum 
electron Lorentz factor to values smaller than "f sa t,e- The 
electrons reaching the most extreme energies will have to 
remain in the acceleration region throughout their life, so 
that they will preferentially cool in the Weibel-generated 
fields, rather than in th e background field B p (at odds 
with the assumption by Kumar et al. ( 2012[ )). By com- 
paring eq. Q and eq. (|14p, we can constrain the combi- 
nation \ c /u vi e B — 0.03! By assuming that the charac- 
teristic transverse wavelength of the Weibel fluctuations 
is A ~ lOc/tjpi (or equivalently, A c / Wp . ~ 10), we can 

estimate €b — 3 x 10~ 3 as the mean value of the mag- 
netic energy fraction in the acceleration region^] This 
is much larger than a, suggesting that the electrons will 
primarily cool due to the self-generated fields (we will 
check this a posteriori) . By balancing the acceleration 
time in eq. ( 14 ) with the synchrotron cooling time in 



the Weibel-generated fields, we infer the cooling-limited 
Lorentz factor in the post-shock framcFj 
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wher e e g = 10 _ e_B._2.5- Remarkably, the limit in 
eq. (26) does not depend on the shock radius. By 
comparing eq. (25) with eq. (26), we see that syn- 



chrotron losses limit the maximum electron Lorentz fac- 
tor ^raax.e = sat,e, 1 sync.e] Only at Small radii, 

11 From now on, we will keep eg as a free parameter. However, 
the combination A 7^ e_g — 0.03 needs to stay fixed, so we expect 
eg to vary by at most a factor of ten, as a result of different choices 



for the transverse scale 3 < A c /^ vi 1^ 30 of the Weibel filaments. 

12 Here, we only consider synchrotron cooling. To include IC 
losses, one would need to generalize the arguments in|Li &: Waxmanl 



losses, one would need to generalize tnc arguments u 
1 2006 1 ; |Piran fc Nakar| |2010[ l; |Li fc Zhao| (2011); 
|2DT2i using tne scalmgs we nave lound m eq. (|l4J 
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whereas for R > 8.5 x 10 17 cm the effect of the flow 



magnetization in eq. ( 25 ) becomes more constraining. 



For electrons, the requirement that the acceleration 
time in eq. ( 14 ) should be shorter than the age of the 
blast wave (~ R/Tc) yields a critical Lorentz factor 



lag 



6.6 x 10 8 < 5 4 4 i?- 7 1/4 



(27) 
than 



which is always larger (so, less constraining) 

lmax,e = ^-^[l sat.ei 1 sync,e] defined above. 

A posteriori, we can check the validity of our initial 
assumption, that the highest energy electrons primar- 
ily cool due to the Weibel-generated field s, rather than 
the ba ckground field B$. As discussed by |Kumar et al.| 
(2012), this is realized if 



£bL b sat > or L ^{^ max ^ e ) 



(28) 



where r L ^{j max ^ e ) _ = (7 mQX , e m e /7 m l ) a 1/2 c/w pi is 
the Larmor radius in the background field Bq, and we 
have taken Lb, sat from eq. (JsJ) as a conservative lower 
limit for the thickness of the region filled with Weibel 
turbulence around the shock (in reality, Lb, sat only ac- 
counts for the extent of the turbulen t re gion on the up- 
stream side of the shock). In eq. (281, the left-hand 
side quantifies the energy lost to synchrotron emission 
as the accelerated particles propagate across a region of 
extent ~ Lb, sat filled with Weibel fields of energy den- 
sity oc es- In contrast, the right-hand side estimates the 
cooling losses due to the background field (with energy 
density oc a) during one Larmor period of the highest en- 
ergy electrons. From the expression for Lb, sat in eq. ([6]), 
we can rewrite the previous inequality as 



£b/c > Jrnax,e/lpk,e 



(29) 



where we have assumed rji n j = 1 in eq. (Km for simplic- 
ity, and we have used that "f p k,e ~ 7o m i/ m e- For our 
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eq. (24 1, and 7 maa;ie from the mini mum between eq. (|25) 



Ipk, 



and eq. (26)), the inequality in eq. (29) is easily satisfied 
for all Lorentz factors < 7 ma s,ej^] 

It is also straightforward to verify that e-BLs.sat ^ 
aR/T during the whole relativistic deceleration p hase , 
a criterion that is even more constraining than eq. ( 28 ) . 
Here, we are comparing the cooling losses in the Weibel 
fields (across a region of thickness ~ Lb. sat) with the 
energy lost to synchrotron radiation in the background 
field, assuming that it permeates the whole GRB shell (of 
thickness ~ R/T, as measured in the post-shock frame). 
This confirms that most of the synchrotron cooling oc- 
curs in the self-generated Weibel fields. 

From the expressions for j S ync,e and j sa t,e derived 
above, we can compute the corresponding synchrotron 
frequencies as a function of the observer's time. We use 
that the observer's time is related to the shock radius 
by t obs = (1 + z)R/2T 2 (e.g., |Waxman||1997| ), where z is 
the redshift. By boosting the emitted comoving frequen- 
cies by a fa ctor of ~ 2T/(1 + z) to obtain the observed 
frequencies (Piran & Nakar 2010), we obtain 
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13 |Kumar et al.| |2012| came to the opposite conclusion, because 
they (incorrectly J employed a value of Lg aat ~ c/wpj that is al- 
most five orders of magnitude smaller than what we find based on 
the results of our PIC simulations. 
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where t i, s = 10 2 t o b s ,2 s - The maximum energy of 
synchrotron afterglow photons will satisfy hv max , e = 



early times (t b s < 4 x 10 s), whereas 



Ghiscllini et al.||20T0~), as 



Acker- 



hl^sync,e at 

hv max , e — hv S at,e at later times. Our results show that 
shock-accelerated electrons cooling in the Weibel fields 
can radiate synchrotron photons up to energies of a few 
GeV in the early phases of GRB afterglows. This ex- 
plains both the GeV photons seen ~ 100 s after the ex- 
plosion in a number of Fermi - LAT GRBs (e 
mann||2010| |De Pasquale|2010 
well as the > 100 MeV photons seen in large numbers up 
to ~ 1000 s after the explosion. 

5.1.2. Stellar Wind 

For the sake of completeness, we also consider the case 
of an adiabatic GRB blast wave propagating in the wind 
of its progenitor star. Let us define A = rriinR 2 , where 
n oc R~ 2 is the particle density in the wind. The Lorentz 
factor of the blast wave will evolve as 



T(R) 



V I6n Ac 2 J 



(32) 



from 

io 15 ^ ,5 4 r ( 7j. 5 ^ 



radius R dec 



6.3 



x 



the deceleration 

~ 1 1 cm up to the non-relativistic 



radius R r , 



11.5 

6.3 x 10 20 ^o,54-4h 1 ! 



cm. Here, we have 



defined A = 3 x 10 n j4n.5 gem^ 1 , as appropriate for the 
dense wind of Wolf-Rayet stars. 
For protons, eq. (16) states that the magnetization of 



the stellar wind constrains the maximum Lorentz factor 
in the wind frame to be 

-1/4, 



up 

' sat.i 



1.3X10 7 E 0M A^ 5 aZl /4 R 



16 



(33) 



where a = 10~ 8 er_8 and R = 10 16 i?i6 cm. Using the 
same argument as in §5.1.1[ acceleration to higher ener- 
gies is possible if the longitudinal coherence length of the 
wind magnetic field satisfies (in the wind frame) 



\coh < 8.4 x 10 13 E^A^a 



-1/2 D l/2 

-8 R 16 Cm 



(34) 



so that the shock can cross a region of length X co h before 
the protons are accelerated up to 7"^ t . For the acceler- 



ation timescale in eq. (14), we have used that the proton 



plasma frequency in the wind is 

w pi ~ 5.7 x 10 4 A{{ 2 5 R^ Hz . (35) 

In the equatorial plane of the stellar wind, the polarity 
of the toroidal magnetic field will alternate with the stel- 
lar rotation period. The half-wavelength of the Parker 
spiral will be Aw ~ 7r(t>w /v ro t)R*, and for typical pa- 
rameters of Wolf-Rayet stars (ratio of wind velocity to 
stellar rotation velocit y vw/vrot ~ 10, and stellar ra- 



dius i?* 
X w ~ 4.4 
eq 



2Oi? , 
10 13 



see 
cm, 



Eichler & Usov 1993), we find 



(34). 



x ±ir ,J cm, which satisfies the condition in 
It follows that acceleration of protons to en- 



ergies larger than in eq. ( 33 1 is generally allowed in the 
equatorial plane of the stellar wind. 

A separate constraint comes from the age of the blast 
wave (~ R/Tc), which limits the maximum proton 
Lorentz factor to be, in the wind frame, 



up 

la 



ge,i 



1.4 x lO^A-fi? 



-3/4 



16 



(36) 
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which is a strict upper limit that does not depend on the 
coherence scale of the circum-burst magnetic field. 

With regards to electrons, the thermal peak of their 
energy spectrum, which we take as a proxy for the low- 
energy end of the electron power-law tail, will be 



ipk,t 



fm, 
2m P 



2.3 xlO^Mfl^ 2 



(37) 



From eq. ( 16 ), the magnetization of the stellar wind con- 



strains the maximum electron Lorentz factor to be 



7sai,r 



9 2 x 10 7 E 1/2 A~ 1/2 a~ 1/4 R' 1 ' 



-1/2 



(38) 



Another co nstr aint comes by balancing the acceleration 
time in eq. ( 14 ) with the synchrotron cooling time in the 



Weibcl fields, which gives 



7s 



(39) 



so that the maximum electron Lorentz factor is set ini- 



tially by ' 



7s 



it increases in radius as oc i? 1 / 3 



(see eq. (39)) up to R ~ 5.0 x 10 17 cm, and for larger 
distances it decreases as oc R- 1 / 2 following 
eq. (38) (but it may be even larger than 7 sa t ie in the 



equatorial plane of the wind, where the field orienta- 
tion chang es wit h the stellar rotation period). As ob- 
served in { 5.1. 1| for the ISM case, the age of the blast 
wave does not play a major role in limiting the maxi- 
mum electron energy. From the expression for j maXte = 
min[7 sa t iE , jsync.e] obtained above, it is easy to verify 
that es/ff> ima x.el lvk.e even in the case of a wind-like 
medium (see eq. (29) for the ISM profile), so the highest 



energy electrons primarily cool in the Weibel-generated 
fields, rather than in the background field. 

From the expressions for 7 S j,„ Cie and 7 sa t,e: we obtain 
the maximum energies of synchrotron photons 



hv -11 F 1/3 A~ 1/e < 
hv sat e ~ 37.1 Ell^e 1 ' 



0,54 c B,-2.5 - 



1/6 (l + Z )- 1 / 3 t- 2 ^GeV(40) 
^ /2 (l + -) 1/2 U 3 A 2 GeV(41) 



The maximum energy of afterglow photons will satisfy 
hv max ,e = hv sync ^ at early times (t obs < 7 x 10 3 s), 
whereas hv max ^ e — hv sat ^ e at later times. Synchrotron 
emission from the external shock of a GRB blast wave 
propagating in a stellar wind is then a valid candidate for 
producing the sub-GeV signature of Fermi-LAT bursts. 

The expression derived above for hv sa t,e can be used 
to interpret the very sharp temporal break observed at 
t obs - 4 x 10 7 s in the light curve of GRB 0607 29, the 
burst with t he latest Chandra X-ray detection. |Grupe| 



et al. (20101 found that the model of a spherical blast 
wave wTEn — Eo j 54 = 1 propagating in a stellar wind with 
^4ii.5 = 0.1 fits the data very well until the temporal 
break at t t, s — 4 x 10 7 s. Following the suggestion by 



Sagi & Nakar (2012), we claim that the temporal break 
at late times (which has been tentatively associated to 
a softening in the 0.3-10 keV spectrum) could be caused 
by the inability of the shock to accelerate X-ray emitting 
particles, due to the limitation given by the wind magne- 
tization. At t }, s ~ 4 x 10 7 s, the shock is still relativistic 
(r ~ 7), so our results are still applicable. By impos- 
ing hv sat ^e — 0.3 keV at t b 3 — 4 x 10 7 s for a burst at 
z = 0.54, as appropriate for GRB 060729, we constrain 
the magnetization of the stellar wind to be a ~ 3 x 10~ 9 . 



5.2. Pulsar Wind Nebulae 

The spectrum of Pulsar Wind Nebulae consists of two 
components, where the low energy component, most 
likely dominated by synchrotron, shows a cutoff at a few 
tens of MeV. The fact that synchrotron emission reaches 
these energies, despite the rapid synchrotron cooling, im- 
plies that particle a cceleration in the nebula i s an ex- 
tremely fast process (de Jager & Harding 19921 . In this 
section, we study the acceleration properties of the ter- 
mination shock in pulsar winds, taking the Crab Nebula 
as our prototypical example. 

Around the equatorial plane of obliquely-rotating pul- 
sars, the wind consists of toroidal stripes of opposite mag- 
netic polarity, separated by current sheets of hot plasma. 
It is still a subject of active research whether the alternat- 
ing stripes will dissipate their energy into particle heat 
ahead of the termination shock, or whether the wind 
remains do minated by Poynting flux till the termina- 
tion sh ock (jLyubarsky fc Kirk||2001| |Kirk fc Skj aeraasen; 
2003). jLyu barsky (2003) showed that, as regards to the 
post-shock properties of the flow, the effect of the stripe 
dissipation at the termination shock is equivalent to their 
annihilation well ahead of the shock. However, when con- 
sidering the synchrotron losses of the accelerated par- 
ticles in the pre-shock fields, it does make a difference 
whether or not the stripes are surviving till the termi- 
nation shock. Here, we only consider the case where the 
stripes are dissipated far ahead of the termination shock. 

The pulsar wind is an electron-positron flow, moving 
towards the termi nation shock wit h a relativistic bulk 
Lyu"barsky||2003[ ) 



Lorentz factor (see 



L 



7o 



sd 



m e c 2 N 



3.7 x 10 A L sd 



38. 5 -^M) 1 



(42) 



where we have assumed complete dissipation of the 
alternating stripes in the wind. Here, L s( i = 3 x 
38.5 ergss 1 is the spin-down luminosity of the 

Crab, and N = 10 40 iV4o s _1 is the total particle flux en- 
tering the nebula, including the long -lived radio-emitting 
electrons (Bucciantini et al. 2011). More precisely, a 



particle flux of 10 4U s _i is a strict lower limit for radio 
emitting electrons, since this estimate only accounts for 
the pairs needed to produce the Crab radio emission at 
frequencies above the ionospheric cutoff - obviously the 
radio spectrum must extend to lower frequencies, there- 
fore the particle distribution goes to lower energies, and 
the inferred number of electrons will increase as the low- 
energy cutoff recedes. On the other hand, one could 
speculate that radio-emitting electrons (whose flux is 
N > 10 40 s~ 4 ) are only injected in the polar flow, and 
that the particle flux in the equatorial region, where the 
electrons emitting in the optical, X-ray and gamma-ray 
bands are ac celerated, will be lower (N — 3 x 10 38 s _1 ), 



as argued by Arons (1998). With these uncertainties in 



mind, we take N = 10 40 s _1 as our reference value. 

Even in the case of complete dissipation of the alter- 
nating fields, the pulsar wind away from the midplane 
will carry a net ordered field, resulting from the fact 
that above and below the equatorial plane the stripes 
are not symmetric (i.e., the widths of the regions with 
opposite fields are not equal). If £ ~ 45° is the inclina- 
tion angle between the rotational and magnetic axis of 
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the Crab pulsar, inferred from fitting the spect rum and 
pulse profil e of the high-energy pulsar emission ( Harding 
et al.|2008|) , and if Aj is the latitude angle measured from 



the midplane, the residu al wind magnetiza tion near the 



gn 

equatorial plane will be (Komissarov|2012 1 



<7(\) 



tan A^ 
tan £ 



(43) 



where we have assumed a <C 1. Based on the results of 
£3.2.2 we require a < 10~ 2 for efficient particle accel- 
eration in electron-positron flows, which corres ponds to 
a < 0.1. T his is th e same constraint found by ISironi & 



Spitkovsky (2011a I with PIC simulations of striped pul 
sar winds, even though they assumed that the magnetic 
stripes survive till the termination shock. In other words, 
we find that the limit on a for efficient Fermi accelera- 
tion does not depend on the location where the magnetic 
stripes dissipate (either in the wind, or at the shock). 

As in the case of GRB afterglows, an important con- 
straint for the maximum energy of shock-accelerated par- 
ticles comes from the flow magnetization. By defining 



a = 10 a_4, eq. (15 1 gives 



Jsat,e( a ) — I' 5 X 10 7 i s d,38.5A r 40 la -4 



1/2 



(44) 

which does not provide any constraint in the midplane 
(where a = 0). Here, we have adopted a local descrip- 
tion of the acceleration process, such that the accelerated 
particles only experience the conditions at one single lati- 
tude. For high-energy electrons, this assumption is likely 
to break, and the acceleration process will be modified 
by the latitudinal dependence of the wind magnetization. 
Such a detailed analysis is beyond the scope of this work. 

Cooling in the self-generated Weibel fields provides an 
additional constraint, independent fro m th e latitude. By 
balancing the acceleration time in eq. ( 13 1 with the syn- 
chrotron cooling time in the Weibel fields, we obtain 



lsync,e — 3.5 X lO S L sd 3g 5 N 4Q 



-1/3-1/3 pl/3 



e B,-2.5^TS,17.5 ■ (45) 



To compute the positron plasma frequency in eq. ( 13 ), we 



have used that the number density ahead of the termina- 
tion shock is nxs = N / (4ir R^ s c) , assuming an isotropic 
particle flux, so that the positron plasma frequency is 



>jj r; 



' 27rnTse 2 
7o m e 



0-11^8.5^0^17.5 Hz (46) 



where the termination shock radius has been written as 

i? TS = 3 X 10 17 i?TS,17.5 cm. 

A third constraint, still independent from the latitude, 
comes from the requirement that the diffusion length of 
the highest energy electrons be smaller than the termi- 
nation shock radius (i.e., a confinement c ons traint). Al- 
ternatively, the acceleration time in eq. ( 13 ) should be 
shorter than i?xs /c, which yields the critical limit 



7eon/, e 



1.9 x 10 7 ^/ 4 38 . 5 7V 40 



1/2 



(47) 



which is generally more constraining than the cooling- 



limited Lorentz factor 7 S 



discussed above. As com- 



pared to the latitude-dependent 7 sa t, e , we see that the 
maximum electron Lorentz factor will be controlled by 
Iconf.e for a_4 < 0.6, and by 7 sa i je for higher latitudes. 



The synchrotron photons emitted by the electrons with 
Lorentz factors 7 sa t, e and 7 co „/. e will have energies 

hv sat , e ~ 96.5 ^ 2 3 8 . 5 A r 4o 2 4 /2 -2.5 i?; ?s,i7.5«-l eV(48) 

^ co „ /ie ~0.17L 2 d! 3 8 . 5 ^ 4 Ve^ /2 _ 2 5 i?^, 17 . 5 keV (49) 

which are apparently too small to explain the X-ray 
spectrum of the Crab, extending to energies beyond a 
few MeV. This conclusion holds even in the extreme 
case 6b ~ lj which might be appropriate in comput- 
ing hv sat e and fw con t e if the stripes survive till the ter- 
mination shock. Moreover, even if we neglect the con- 
straint given by fw con f, e , we find that X-ray photons with 
hv > 1 keV can be produced only within an equatorial 
wedge of angular extent A^ < 10~ 5 (assuming an incli- 
nation £ ~ 45°). The energy flux in this sector is a tiny 
fraction ~ 2a ~ 2 x 10~ 5 of the total energy flow in the 
wind (assuming the energy flow follows the Poynting flux 
of a split monopole cx cos 2 Af), which is insufficient to 
power the Crab X-ray emissi on, that consumes ~ 10% o f 
the pulsar spin-down power ( jAtoyan &: Aharonian|19"96 ). 
We conclude that Fermi acceleration at the termination 
shock of PWNe is not a likely candidate for producing 
X-ray photons via the synchrotron process, and valid al- 
ternatives should be investigated. 

One possibility - magnetic dissipation of the striped 
pulsar wind in and around the shock front itself - has 
been extensively studied, with the conclusion that parti- 
cle acceleration along extended X-lines formed by tearing 
of the current sheets may contribute to the flat particle 
distribution (with spectral index p ~ 1.5) required to ex- 
plain the far infrared and radio spectra of P WNe (e.g., 



Lyubarsky 2003 Sironi & Spitkovsky 2011a). However 



further acceleration to gamma-ray emitting energies by 
the Fermi process cannot occur in the transverse shock 
that terminates the pulsar wind, if particle scattering de- 
pends only on the self-generated turbulence. 

Yet, the steady-state hard X-ray and gamma-ray spec- 
tra of PWNe do look like the consequences of Fermi accel- 
eration - particle distributions with p ~ 2.4 are implied 
by the observations. In this regard, we argue that the 
wind termination shock might form in a macroscopically 
turbulent medium, with the outer scale of the turbu- 
lence driven by the large-scale shear flows in the nebula 

omissarov &: Ly ubarsky||2004| |Del Zanna et al. 1120041 
|Camus et al.||2009| ). If the large-scale motions observec 
in MUD simulations of PWNe drive a turbulent cascade 
to shorter wavelengths, back-scattering of the particles in 
this downstream turbulence, along with upstream reflec- 
tion by the transverse magnetic field of the wind, might 
sustain Fermi acceleration to higher energies. Yet, the 
turbulent cascade has to work in such a way as not to 
disturb the polar ization of the nebula, which looks rather 
cleanly toroidal ( |Wilson|[l972"l |Schmidt et al.|1979 |. 

Another "external" influence ot reconnection on the 
shock structure, that might lead to particle acceleration 
to higher energies, may be connected to the accelera- 
tor behind the re cently discov ered gamma-ray flares in 
the Crab Nebula ( Abdo 2011). If the stripes decay well 
ahead of the shock, the wind has a "Mexican hat" mag- 
netic geometry, with oppositely wound toroidal magnetic 
field in the northern and southern hemispheres, sepa- 
rated by the equatorial current sheet. Tearing of that 
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current sheet can create radial spokes of current, with 
radially extended X-lines. Runaway acceleration of elec- 
trons and positrons at those X-lines, a linear accelerator, 
injects energetic beams into the shock, with the mean 
energy per particle approaching the whole op en field line 
voltage, > 10 16 V in the Crab ( |Arons||2012j ). This high- 
energy population can drive cyclotron turbulence when 
gyrating in the shock-compressed fields, and resonant ab- 
sorption of the cyclotron harmonics can accelerate the 
electron-positron pairs in a broad spectrum, with max- 
imum energ y again comparable to the whole open field 
line voltage (Amato & Arons 2006). 

6. SUMMARY AND DISCUSSION 

In this work, we have investigated by means of 2D and 
3D PIC simulations the physics and acceleration proper- 
ties of relativistic perpendicular shocks, that propagate 
in electron-positron or electron-ion plasmas with moder- 
ate magnetizations (0 < a < 10 _1 ). As a function of the 
magnetization a and the pre-shock Lorentz factor 70, we 
have explored the efficiency of the Fermi process (i.e., the 
amount of particles and energy stored in the non-thermal 
tail) and the rate of particle acceleration. 

We find that the Fermi process is suppressed in 
strongly magnetized perpendicular shocks, in agreement 
with the conclusions by SS09 an d SS11 and with ear- 
lier ID simulations perfo rmed by Langdon et al. ( 1988 1 
and Gallant et al. (1992). Due to the lack ot sufficient 
self-generated turbulence, the charged particles are con- 
strained to move along the field lines, which are advected 
downstream from the shock. The Fermi process, which 
requires repeated crossings of the shock, is then inhibited. 
As a result, the post-shock particle spectrum is consis- 
tent with a Maxwellian distribution, for magnetizations 
a > 3 x 10~ 3 in electron-positron flows and a > 10 -4 in 
electron-ion flows. We have tested that the lack of par- 
ticle acceleration in high-cr shocks is a solid result, con- 
firmed by means of 3D simulations and validated with 
different choices of the mass ratio (up to rrii/m e — 1600, 
close to the realistic value mi/m e ~ 1836). 

Weakly magnetized flows (a < 10~ 3 in electron- 
positron flows; (T < 3 x 10~ 5 in electron-ion plasmas) 
are mediated by the Weibel instability, as in unmagne- 



positrons, equal to the electron plasma frequency uj pc ) 
and £ m ax/lo m iC 2 — 0.25(w p i<) 1 I 2 in electron-ion flows 
(here, w p ; = \J m e /rrii w po is the ion plasma frequency). 
The relation e max oc i 1 / 2 , rather than the commonly as- 
sumed Bohm scaling e max oc t, is a natural conseq uence 



tized shocks (|Spitkovsky||2008a|b| |Martins et al.||2009 
Haugb0lle|2Olip . The instability, seeded by the counter- 
streaming between the incoming flow and a beam of 
shock-reflected particles that propagate back into the 
upstream, produces strong small-scale magnetic turbu- 
lence that governs the Fermi process. As a result, the 
post-shock spectrum in weakly magnetized shocks shows 
a prominent non-thermal tail of shock-accelerated parti- 
cles. Regardless of a (in the regime of weakly magnetized 
flows) or 70 (provided that 70 > 10), the non-thermal tail 
contains ~ 1% of particles and ~ 10% of flow energy, and 
its power-law slope is p ~ 2.5. In electron-ion shocks, the 
two species enter the shock nearly in equipartition, as a 
result of strong energy exchange ahead of the shock me- 
diated by the time-varying Weibel fields. As a result, 
electrons and ions have similar acceleration efficiencies. 

For weakly magnetized shocks, the non-thermal 
tail of shock-accelerated particles stretches in time 
to higher and higher energies, following the scal- 
ing emax/lo^iC 2 — 0.5(w p it) 1 / 2 in electron-positron 
plasmas (here, w p i indicates the plasma frequency of 



of the small-scale nature of the Weib el turbulence (Kirk 
fc Reville|20l0| |Plotnikov et al.|2012"| . The coefficienTof 
proportionality (i.e., ~ 0.5 for electron-positron, ~ 0.25 
for electron-ion) provides an estimate of the combina- 
tion esA c / Wpi , as shown in eq. (Ej). Here, es is the frac- 
tion of flow energy converted into magnetic fields in the 
shock region, and A c / Wp . is the transverse scale of the 
Weibel modes in units of the plasma frequency c/w p j. We 
find X c f Upi ~ 10, which implies es — 10 -2 in electron- 
positron flows and cb — 3 x 10~ 3 in electron-ion plasmas. 
Our estimates for es do not depend on a (in the regime of 
weakly magnetized flows) or 70 (provided that 70 > 10). 

By measuring the time t(e) needed to accelerate the 
particles up to a given energy e, our simulations are im- 
plicitly providing an estimate for the diffusion coefficient 
D ~ c 2 t in relativistic perpendicular shocks. We find 
that D ~ 4 c c/uj p i(e/jomiC 2 ) 2 in electron-positron flows, 
and D ~ 16 c c/w p i(e/7oTOiC 2 ) 2 in electron-ion plas- 
mas, regardless of the flow magnetization a or the bulk 
Lorentz factor 70. This clearly differs from the so-called 
(and widely used) Bohm scaling, where Db ~ cr^, a {£) ~ 
er~ 1//2 c c/o; pi (£/7omjC 2 ). The two estimates differ by a 
factor of D/D B — 4 a 1 ' 2 (e/7o m^c 2 ) in electron-positron 
flows and D/Db — 16 a 1 ^ 2 (e/^fo m iC 2 ) in electron-ion 
plasmas. For small a, this shows that the Fermi process 
in relativistic shocks may be even faster than Bohm, if 
we calibrate the Bohm scaling with the pre- shock fieldj^ 

The scaling e max oc t 1 / 2 breaks down at an energy 
e sa t that depends on the flow magnetization a. The 
thickness of the upstream region populated with Weibel 
filaments scales as the Larmor radius of the particles 
at the low-energy end of the non-thermal tail, giving 
Lb, sat ~ c -1 / 2 c/w p i. When the diffusion length of the 
highest energy particles (~ D/c) reaches Lb, sat, further 
evolution to higher energies is inhibited, due to the lack of 
sufficient magnetic turbulence to oppose advection down- 
stream with the ordered field. At this point, the shock 
reaches a steady state, and the particle spectrum satu- 
rates with an upper energy cutoff £ S at/lo m iC 2 — 4er~ 1//4 
in electron-positron shocks and e S at/7o TO iC 2 — 2 cr -1 / 4 in 
electron-ion flows. The saturation in the particle spec- 
trum has been confirmed with 3D simulations. Accel- 
eration to higher energies (beyond £ sa t) might be gov- 
erned by pre-existing upst ream turbulence, as argued by 
Sironi & Goodman (2007), provided that sufficient tur- 
bulent power exists on scales ~ L B , S at — 10 13 fjlg^ 2 cm 
(as measured in the post-shock frame). 

Our results can provide physically-grounded inputs for 
models of non-thermal emission from a variety of astro- 
physical non-thermal sources, with particular relevance 
to GRB afterglows and t he te rmination shock of PWNe. 
As we have described in §5.1[ we find that the accelera- 
tion of protons by the external shocks of GRB afterglows 

14 When the Bohm scaling is computed with the self-generated 
fields, in the formulae above we should replace cr with 2eg . 
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is too slow (given that t oc e 2 ) to explain the extreme 
Lorentz factors of UHECRs. In the early phases of GRB 
afterglows, electrons can be accelerated up to Lorentz 
factors ~f sy nc,e ~ 10 7 before suffering catastrophic syn- 
chrotron losses in the Weibel-generated fields (as opposed 
to the backg round pre-shock field, as argued by [Kumar] 
et al. (2012)). Their synchrotron radiation can produce 
the ~ (jev photons detected by the Fermi telescope 



e.g., AckermannpOlO De Pasquale 2010] Ghiscllini et al. 

TTT 

J It 



2010 1 in early GRB afterglows. A more comprehensive 



comparison with the obs ervations, including t he effect 
of IC losses (as done by Li fc Waxman|[2006| |Piran fc 



20TTT 



NakarpOTO] |Barniol Duran & Kumar||201H |Li fe Zhao 
|L 



emome 



2012] |Sagi & Nakar||2012p 5ut without a 



self-consistent model for the Fermi process in GRB af- 
terglows), is deferred to a future study. 

As regards to PWNe, Fermi acceleration in the mag- 
netized transverse relativistic shock that terminates the 
pulsar wind has been believed, for more than 40 years, 
to be the mechanism behind the conversion of the wind 
energy into the non-thermal particle spectra inferred in 
PWNe. In 95.2] we have shown that Fermi acceleration 
in the turbulence self-generated at the termination shock 
is not able to accelerate the particles to sufficiently high 
energies to explain the X-ray and gamma-ray emission 
from the Crab Nebula. To explain the observed spectral 



signatures of PWNe, we suggest that the turbulence re- 
quired for the Fermi process might be generated either 
in the nebula, by large-scale motions dCamus et al.|200"9| 
casca ding down to smaller wavelengths ( Bucciantin i et al. 
201 1[), or at the shock itself, by cyclotron waves ( jAmato 
& Arons 2006] ) excited by a relativistic pair beam that 
was pre-acc eTerated in t he pulsar wind via magnetic re- 
connection ( Arons|[2012 1. 
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APPENDIX 
DEPENDENCE ON THE MASS RATIO 



In Fig. 14 we analyze the dependence of the post-shock spectrum in electron-ion shocks on the mass ratio mi/m e , 
up to mi/m e — 1600, that approaches the realistic value mi/m e ~ 1836. The left column shows the results for 
a = 3 x 10~ 4 , and the right column refers to a = 10~ 5 . We find that the shape of the post-shock spectrum and the 
temporal scaling of the maximum energy of shock-accelerated particles are nearly insensitive to the mass ratio m,/m e . 

For u = 3x 10~ 4 (left column in Fig. 14), the spectrum of ions (top panel) and of el ectr ons (bottom panel) shows 
no evidence for particle acceleration, provided that mj/m e > 100. As discussed in [4.2 the electrons populate a 
Maxwellian distribution, whereas the ion spectrum is not entirely consistent with a Maxwellian, due to incomplete ion 
thermalization at the shock front. As the mass ratio decreases (blue line for rrii/m e = 25, purple line for rrii/m e — 6.25), 
the peak of the electron Maxwellian systematically shifts to higher energies, suggesting an increase in the efficiency of 
ion-to-electron energy transfer ahead of the shock, for smaller m,i/m e . Yet, the location of the thermal peak in the 
electron spectrum only varies by a factor of ~ 3, as the mass ratio changes from mi/m e — 6.25 up to mi/m e — 400. 
Efficient electron heating ahead of the shock requires strong Weibel turbulence, which can be excited if the upstream 
electrons are unmagnetized. This requires that the initial electron magnetization cr e> o = (mi/m e )a should not be too 
large, a constraint that, for fixed u, can be satisfied more easily for smaller mi/m e . This explains the trend in the 



electron thermal peak seen in the bottom panel of Fig. 14 (left column) 



The increase in the number of non-thermal electrons as the mass ratio decreases (bottom left panel of Fig. 14 ) is 



then a consequence of the heating experienced by the incoming electrons in low-mi/m e shocks. If the electrons are 
entering the shock with a larger energy, they are more likely to be injected into the acceleration proce ss, w hich explains 
the trend in the normalization of the electron non-thermal tails seen in the bottom left panel of Fig. [14] Yet, even for 
mi/m e — 6.25, namely the case that shows the largest acceleration efficiency, the maximum electron Lorentz factor 
stops growing after uj p [t ~ 600 (purple line in the inset (b) of the bottom left panel). 



For a = 10 -5 (right column in Fig. 14 1, the post-shock spectrum of both ions (top panel) and electrons (bottom 
pan el) shows a prominent non-thermal tail, extending in time to higher and higher energies. In the right column of 
Fig. 14 we plot the spectra at uj pl t = 562 with solid lines, and at Lo pl t = 1890 with dashed lines. For both times, the ion 



and electron spectra are nearly identical across a wide range of mass ratios: at relatively early times (u! pl t = 562), we 
can follow the physics of electron-ion shocks with large mass ratios (mj/ro e = 100 in green, mi/m e = 400 in yellow, and 
the nearly realistic case rrii/m e = 1600 in red); when evolving the shock to longer times (w p ji = 1890), we are limited 
to relatively small mass ratios (rrii/m e = 6.25 in purple, mi/m e — 25 in blue, mi/m e = 100 in green). The shape 
of the spectra (for both ions and electrons) is remarkably similar across the whole range of mass ratios we explore, 
with the only difference being a moderately larger degree of electron heating for smaller mi/m e , as explained above 
for a — 3 x 10~ 4 shocks. Yet, the difference in the electron peak energy between rrii/m e = 6.25 and rrii/m e = 100 at 
ujpit = 1890 is less than a factor of two, and it is even smaller between rrii/m e = 100 and mi/m e = 1600 at oj pl t = 562 



As time progresses, the non-thermal tail of both ions (top panel in the right column of Fig. 14) and electrons (bottom 
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Fig. 14. — Comparison of downstream particle spectra (upper panel for ions, lower panel for electrons) between different mass ratios, 
from m,i/m e = 6.25 up to m;/m e = 1600. The left column refers to a = 3 X 10~ 4 , and the right column to a = 10 -5 . In the insets, we 
plot the temporal evolution of the maximum particle energy. 

panel) systematically extends to higher energies, as exemplified by the case with mi/m e = 100 (compare the solid 
green lines for oj v \t = 562 to the dashed green lines for w v it = 1890). The tem poral evolution of the maximum ion and 
electron Lorentz factor is plotted in the insets of the right column of Fig. 14 (inset (a) for ions and (b) for electrons), 



showing that the scaling as oc (wpji) 1 / 2 is realized at late times independently of the mass ratio. Based on this evidence, 
we argue that the results presented in the main body of the paper, both in terms of acceleration efficiency and of 
energization rate, are going to hold even for the realistic mass ratio mi/m e ~ 1836. 

DEPENDENCE ON THE FIELD ORIENTATION RELATIVE TO THE SIMULATION PLANE 

In the main body of the paper, we have preferentially employed 2D computational boxes with the pre-shock magnetic 
field Bq oriented perpendicular to the simulation plane (i.e., out-of-plane fields). Here, we demonstrate that this 
configuration is in excellent agreement with 3D simulations. 

In Fig.[l5j we show the post-shock spectrum for different geometries of the upstream background magnetic field. We 
compare 2D simulations with in-plane fields (green) , 2D simulations with out-of-plane fields (black) , and 3D simulations 
(red), for both electron-positron flows with a = 10~ 3 (left column) and electron- ion flows with mj/m e = 25 and 
a = 3 x 10~ 5 (right column). The disagreement between 3D simulations and 2D runs with out-of-plane fields at low 
energies (below the thermal peak) is a mere consequence of the different adiabatic index, which changes the analytic 
formula of a Maxwellian distribution. In fact, for 2D simulations with out-of-plane fields, the plasma is constrained 
to move in the plane orthogonal to the field, which coincides with the simulation plane. In contrast, in 2D runs with 
in-plane fields, these two planes are not degenerate, the particle velocities can sample all three dimensions, and the 
adiabatic index is the same as in 3D (which explains the agreement at low energies between the green and red lines). 

At high energies (i.e., in the non-thermal tail), 2D simulations with out-of-plane fields are in good agreement with 
3D simulations, whereas 2D runs with in-plane fields systematically underestimate the acceleration efficiency. 3D 
simulations and 2D runs with out-of-plane fields yield consistent results as regards to both the acc eler ation efficiency 
and the temporal evolution of the maximum energy of accelerated particles (see the insets in Fig. 15). In summary, 
the 2D simulations with out-of-plane fields employed in the main body of the paper provide an excellent description 
of the 3D acceleration physics of relativistic perpendicular shocks. 
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Fig. 15. — Comparison of downstream particle spectra among different geometries of the upstream background magnetic field. We 
compare 2D simulations with in-plane fields (green), 2D simulations with out-of- plane fields (black; this is the configuration adopted in the 
main body of the paper), and 3D simulations (red). The left column refers to an electron-positron shock with magnetization a = 10 — 3 , 
whereas the right column shows the ion and electron spectra of an electron-ion shock (mi/m e = 25) with a = 3 X 10 -5 . In the insets, we 
plot the temporal evolution of the maximum particle energy. 

The artificial s uppression of the a cceleration efficiency in 2D simulations with in-plane fields may be explained by 



Jones' theorem (Jones et al. 19981, that charged particles cannot move farther than one Larmor radius from the 
plane defined by the background magnetic field (along y) and the z direction orthogonal to the simulation domain. 
Equivalently, we could argue that the diffusion of particles back into the upstream, that is required for injection into 
the Fermi process, is likely to be suppressed in the case of in-plane fields, since the growth of the Weibel magnetic 
fields (preferentially along z, for 2D runs) cannot efficiently compete with advection by the ordered background field 
(directed along y). In contrast, if the pre-shock field is initialized orthogonal to the simulation plane, its effects can 
be counteracted more easily by the Weibel modes, and injection into the Fermi process would be facilitated. 

Finally, we remark that the difference seen in Fig. |T5| between the acceleration capabilities of 2D shocks with in-plane 
and out-of-plane fields is most dramatic for magnetizations close to the tr ans ition between poorly-accelerating and 



efficiently-accelerating shocks (this is what motivated the choice of a in Fig. 15 1. For magnetizations that are a factor 
of ~ 3 larger (or smaller) than this critical boundary, the post-shock spectra m 2D simulations are essentially the same 
regardless of the field orientation relative to the simulation plane. 
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